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ON MELTING AND FREEZING FOR THE 2D RADIAL STEFAN 

PROBLEM 

MAHIR HADZIC AND PIERRE RAPHAEL 


Abstract. We consider the two dimensional free boundary Stefan problem de¬ 
scribing the evolution of a spherically symmetric ice ball {r < A(f)}. We revisit 
the pioneering analysis of 1201 and prove the existence in the radial class of finite 
time melting regimes 


A(£) 


(T - t )l/2 e -TT-'/|log(T-t)|+0(l) 
(c + o(i)) (T ~ t)j T +1 , fce N* 

|log(T-t)|“2?T 


as t ■ 


T 


which respectively correspond to the fundamental stable melting rate, and a 
sequence of codimension k £ N* excited regimes. Our analysis fully revisits a 
related construction for the harmonic heat flow in l42l by introducing a new and 
canonical functional framework for the study of type II (i.e. non self similar) blow 
up. We also show a deep duality between the construction of the melting regimes 
and the derivation of a discrete sequence of global-in-time freezing regimes 

Aoo — A(f) ~ | 1 : 'j i- c m* as t +°° 

( t fc (logt ) 2 ’ h ^ 


which correspond respectively to the fundamental stable freezing rate, and ex¬ 
cited regimes which are codimension k stable. 


1. Introduction 

1.1. Setting of the problem. We consider the classical two dimensional one-phase 
Stefan problem on an external domain. The unknowns are the moving domain 
12(f) C K 2 and the temperature function u : 12(f) — >• M which evolve according to: 

{ dtu — Au = 0 in 12(f) 

d n u = V d n(t) on <9f2(f) (1.1) 

u = 0 on <9f2(f) 

where Vqq^ stands for the normal velocity of the moving boundary <912 (£0 The 
temperature u may either be assumed to be positive initially in 12(0), in which 
case the maximum principle and the Dirichlet boundary condition ensure that it 
will remain positive in 12(f), or on the contrary the data may be undercooled with 
initially non positive temperature in some regions in space. The cavity represents 
a circular block of ice kept at constant temperature u = 0. If the cavity vanishes at 
a later time we refer to this process as melting and if it expands, we refer to it as 
freezing. 


2 For any given parametrization y(f, •) : R —^► R 2 of <9f2(f), the normal velocity is given by the 
formula Vgn(t) = dt'y • n, where n stands for the outward pointing unit normal with respect to 
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1.2. Cauchy theory. There exists a large number of works pertaining to the ques¬ 
tions of existence, uniqueness, regularity, and well-posedness for the classical Stefan 
problem. Weak solutions were first defined and shown to exist in [ 22] and their prop¬ 
erties were further studied in many works, see [m m Euini mi ss] and references 
therein. The classical Stefan problem lends itself to a different notion of a weak 
solution, the so-called viscosity solutions. For an overview of the seminal works on 
the regularity theory for such solutions, we point the reader to [3] and references 
therein, and for the basic existence and uniqueness theory to [231 El]. In the class 
of classical solutions, first existence results are traced back to |18l [30 ]. For well- 
posedness results in energy-based Sobolev-type spaces see mm and in L p -type 
spaces pEi [39]. 

In this article, we work in a radially symmetric situation and therefore the Cauchy 
theory is simpler. The domain Q(t) is given by 

Q(t) = {x € M 2 ; |x| > A(t)}, 


and the Cauchy problem (II.Ill becomes: 

Ut — u rr — A u r = 0 in fl(f) 
< u r (t, X(t)) = -X(t) 
u(t, A (t)) = 0 
„ u( o, ■) = Uo, A(0) = A 0 . 


It is well posed in H 2 : for all (tto, Ao) € H 2 x R*j_ with uq radially symmetric, there 
exists a unique solution (u(t), A(t)) € C([0, T), H 2 (Q)) x C 1 ([0, T), M!j_) to (II.2[) . and 


T < Too implies 


(jim \\ u (t)\\ H *m)) 


Too or 


lim Aft) = 0 

t->T 


We recall the classical proo0 in Appendix [Dj A simple integration by parts using the 
boundary conditions (II. 2D . see (ID.4p . ensures the uniform control of the Dirichlet 
energy: 


'M>A (t) 


\\7u(t,x)\ 2 dx < 


'|x|>A(0) 


\\7uo(x)\ 2 dx 


(1.3) 


in the melting regime A < 0. 


1.3. Main results. The main result of this paper is the existence of a discrete 
sequence of melting rates with k nonlinear instability directions, k S N. 

Theorem 1.1. [Melting dynamics] There exists a set of data (uo,Ao) in H 2 x 
Ml, with uq arbitrarily small in H 1 , such that the corresponding solution (tt, A) € 
C([0,T), H 2 ) x C 1 ([0, T), M_|_) to the exterior Stefan problem H1.2\) melts in finite 
time 0 < T = T(uo, lo ) < oo with the following asymptotics: 

1. Stable regime: the fundamental rate is given by 

A(t) = (T - t) x / 2 e~ ^\/l log ( T -*)ITCW(i) (1.4) 

and is stable by small radial perturbations in H 2 ; this regime corresponds to positive 
data uo > 0. 

2. Excited regimes: the excited melting rates are given by 

k +1 

A(i) = (c*(« 0 , Ao) + ot-»r(l)) (T ~* )a m , k € N*, (1.5) 

|log(T - t)| 2 fc 

2 see for instance 00 and references therein for a Cauchy theory in the class of Holder spaces. 
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for some c*(ito, Ao) > 0; it corresponds to undercooled initial data lying on a locally 
Lipchitz H 2 manifold of codimension k. 

3. Non concentration of energy: In all cases, there exists u* € H 1 such that 

lim||Vu X { N >A W} -Vu*|| L2(R2) = 0. (1.6) 

Comments on the result. 

1. Previous works. Let us recall that a melting scenario with asymptotics on the 
melting rate is first proposed in the pioneering work [ 20] @. The analysis involves 
first a change of variables similar to [21] , the meaning of which for the Stefan prob¬ 
lem is not obvious, and a complicated matching procedure. In a series of recent 
works on both energy critical and supercritical parabolic and dispersive problems 
[331 hd sa H M H, a different angle of attack for the construction of type II 
blow up bubbles is proposed. The strategy consists of two steps: construction of a 
high order approximate solution based on the expansion of the blow up profile with 
respect to a well chosen small parameter, and high order Sobolev energy estimates 
which allow for sufficient decay to close the estimates. The strength of this approach 
lies in its robustness as it in particular applies to Schrodinger and wave problems 
which are more delicate to handle. Its weakness however is that it requires the 
control of a large number of derivatives to get sufficient decay which may not be 
sharp in the parabolic setting. 

2. Role of the dimension. This paper deals with the case d = 2, which is the energy 
critical case, but the higher dimensions d > 3 could be treated by an entirely anal¬ 
ogous approach. In fact, the case d = 2 is the most complicated case, displaying 
small logarithmic gains only and strong coupling between the various components 
of the solution^. 

The main novelty of the proof of Theorem 11.11 is the derivation of a sharp func¬ 
tional setting for the construction of type II, i.e. non self similar blow up bubbles, 
here applied to a melting problem. Our analysis relies on new weighted energy 
bounds with a degenerate Gaussian weight based on the spectral decomposition 
of the leading order linear operator after renormalization. This is conceptually a 
continuation of the Giga-Kohn approach m to the type I blow up in the energy 
subcritical range. Our new set of estimates simplifies both the derivation of the 
approximate solution and the closure of the nonlinear energy bounds by using in an 
optimal way the dissipative structure of the problem, see in particular (|3.48p . The 
existence of a degenerate resonance leading the blow up rate is reminiscent of the 
derivation of the celebrated "log-log law" for the mass critical nonlinear Schrodinger 
equation, see [23 ES E21 El] • A recent series of works by Merle and Zaag [351136] [33 
suggest that this approach may be of great interest for dispersive wave like problems 
as well. 

Moreover, our approach applies as well to the construction of solutions that asym- 
potically converge to the solitary wave (it = 0, A = const > 0), here meaning freez¬ 
ing, and the proof exhibits a deep duality for the derivation of the melting/freezing 
rates. 

^A different version of (11.411 is computed in 1261 . but a correction is suggested in m- A similar 
problem occurs in mi, see the correct law in l44l . 

4 see for example (13.321) . 
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Theorem 1.2. [Freezing dynamics] There exists a set of data (uq, Ao) in H 2 x R+, 
arbitrarily sm.all in H 1 . such that the corresponding solution (u, A) € C([0, T), Ft 2 ) x 
C 1 ([0, T),R + ) to the exterior Stefan problem \1.2\) exists globally-in-time. Further¬ 
more, let 



(1.7) 


then the solution exhibits an asymptotic freezing 

lim A (t) = Aqo > 0 




with the following asymptotics: 

1. Stable regime: the fundamental freezing rate is given by 



c(uq, Aq)(1 + Ot_>+oo(l)) 


( 1 . 8 ) 


for some c(uq, Ao) > 0; it is stable with respect to small well localized smooth radial 
perturbations; this regime contains negative data uq < 0. 

2. Excited regimes: excited melting rates are given by 



for some Cfc(uo,Ao) > 0 and correspond to superheated well localized initial data 


lying on a locally Lipchitz manifold of codimension k in some well localized norm. 

3. Energy asymptotics: in all cases, the Dirichlet energy dissipates at the rate 


dk{u 0) A 0 )(l T Of,— >.-|-oo(l)) 

t k+1 log t 


l|Vw(i)|| L 2 (n(t)} 


( 1 . 10 ) 


k e N. 


for some c4(uo, Ao) > 0. 

Comments on the result. 

1. More melting regimes. In the setting of the Stefan and Hele-Shaw problems, the 
authors consider in m a melting scenario for the one-phase Stefan problem out¬ 
side a fixed domain containing the origin and kept at a pre-fixed non-negative and 
nontrivial temperature, acting as an effective heat source. The liquid thus expands 
to infinity for positive initial data and an asymptotic rate of expansion for the free¬ 
boundary radius is obtained. Note that this situation is quite different from our 
setting as there is no such heat source in our case, and the freezing/melting process 
is driven entirely by the choice of initial conditions. 

2. Solitary wave regimes. A non trivial global-in-time dynamics with convergence to 
the solitary wave similar to that described by theorem 11.21 has been derived in other 
critical settings, see for example m, Ea, m ■ The quantized convergence rates 
with logarithmics corrections are reminiscent of some classical nonlinear dynamical 
systems scenarios, see for example m- 

The first main open problem following this work is the understanding of the full 
non radial stability of the free boundary problem in the stable melting regime k = 0 
which should be amenable to our approach. Let us mention that a related problem 
in the context of evaporating drops was recently studied and solved in the setting 
of a self-similar collapse in the very nice work [7]. The second main open problem 
is to give a complete description of the flow for small initial data, and here the 
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constructions and underlying functional framework of Theorem 11.11 and Theorem 
H72l will be essential steps. 

Acknowledgements. P.R is supported by the Institut Universitaire de France and 
the ERC-2014-CoG 646650 SingWave. M.H. kindly acknowledges the hospitality of 
the Laboratoire J.A. Dieudonne, Universite de Nice Sophia-Antipolis, where part of 
this work has been carried out. 


Notations. We denote the ball of radius K by 

B K (R d ) := {x € R d ,\x\ < K} 

and set 

A := yd y . 

For any a > 0 we denote the external domain: 

fl a := {x € M 2 | |x| > a}. 

When a = 1, we shall simply write fii = fi. We define the weight 


P±{z) ■= e ± 2 


and the scalar product on by 


poo 

(f,g)b,± = / 

J Vb 


fgp±zdz 


and the associated norms 


u L\ + = 

p,b,± 


'Vb 


u~p±zdz) , Hull^i =( 11 ^ 2 ^ 112,2 + \\u 


p,6,± 


p,b ,:t 


I L 2 


p,6,± 


and 


\ u \\H 2 ,, 

p,o,± 


[2 + ||9 z u||?2 + 11 ^ 11 7 2 

P, b,± " L p,b,± L p,b,± 


We define for b > 0 the Hilbert space 

Hp b ,± = {u '■ Qy/b —> K, u radial with |M|#i ft± < +oo and u(Vb ) = 0} (1.11) 

equipped with the scalar product (■, -)k,±, and for 6 = 0: 

Hp 0i= t = {u : M 2 — > R, u radial with ||rt||jp () ± < +oo} 

equipped with the scalar product (•, -)o. Similarily, we define the renormalized 
quantities: 

f°° . b\y \ 2 

(/, g)b,± = j fgpb,±ydy , p b ,± = e* 2 

and the norms 

Uplift,± = (J^ v 2 Pb,±ydy^j , |M|#i ± = IMIm=+||<9jHIm=, \\v\\ H 2 ± = IM|tfi ± + ||Au|| b)± , 

and the Hilbert space 

Hi + = {v : H —>• M, v radial with ||u||^i < +oo and u(l) = 0}. 

o,=b 


( 1 . 12 ) 


^Observe that u(Vb) = 0 ensures that b can be isometrically embedded into Hp 0 by consid¬ 
ering the map u ul z>y /j;. 
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We define the sequence of numbers: 


«0 := 0 , 


a,- : = 


3 -1 

E 

i =0 


3 ~ i 


for j > 1. 


(1.13) 


Throughout the paper, summations over 0 < j < k — 1 are empty for k = 0. 


1.4. Strategy of the proof. Problem CE2D is invariant under an energy critical 
scaling: if (u,X) solves (| 1.2 [) . then so does 

:= A p (t) := —, (1.14) 

[1 

and the scaling leaves the Dirichlet energjj^ unchanged. We therefore renormalize 
the flow 

ds 1 x . 

«(t.O = »(»,!/), S = ^y. y ~ (L15 > 

and obtain the renormalized equation with a /heed boundary: 

{ d s v — — An = 0, y > 1: 

v(a, 1) = 0 (1.16) 

U ?/M) = -T 

step 1 Perturbative spectral analysis. We start with the description of melting 
regimes. Let in a first approximation 

b = — 0 < b <C 1, 

A 

then the linear operator driving (ll.lbp is 

Tib = —A + 6A, u(l) = 0. 

Our first new input is to diagonalize this operator in a suitable Hilbert space. Indeed, 

b\y \ 2 

Tib is self adjoint with respect to the measure e 2 dy, and has up to a shift 
compact resolvent in the Hilbert space and hence discrete spectrum. However, 

the limit b —>■ 0 is singular in the sense that the limiting operator is the Laplacian 
with resonant eigenmodes and continuous spectrum. After renormalization, we may 
equivalently consider 

Hb = —A + A, v{Vb) = 0 

which formally is a deformation of the standard harmonic oscillator, but with a 
singular boundary condition. We claim that in these renormalized variables, the 
operator Hb can be diagonalized using a perturbative Lyapounov-Schnridt type ar¬ 
gument in the weighted Hilbert space b _. The first k eigenvalues are given for 
0 < b < b*{k) by 

X btk = 2k+ ij2_ + o 

with a corresponding asymptotic expansion of the eigenmodes, see Proposition 12.31 
After renormalization, this produces a family of eigenvectors 

Tdbdb,k = bXb t kr]b,k, Vb,k{ 1) = 0. 

step 2 Approximate solution and modulation equations. We now freeze 

A s 

-— = a, 0<a<Sl 

A 

®which is dissipated from d in the melting regime. 


(1.17) 
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and rewrite the renormalized flow 


d s v + T~LbV + (a — b)Av = 0 , v(s, 1) = 0 , d y v(s, 1 ) = a 


for a parameter b which will be chosen later. We look for an approximate solution 
of the form 

k 

Qp(s,y) =^2b j (s)ri b{s ) !j (y). 

3 =0 


After projection onto each eigenmode, we obtain the leading order dynamical sys¬ 
tem: 


(bj) s + bbj\h t k + ( | i og 6 | J + TT^ — 0 ) 0<j<k 
$ = b s + 2 b(a — b ). 


(1.18) 


This sytem is complemented by the renormalized nonlinear free boundary conditon 
d y v(l) = a which forces the leading order relationship 


a = 




+ lot 


with (aij)o<j<k given by (11.131) . It remains to chose b(a ) which is done through the 
choice 

$ = 0 (1.19) 

which will be motivated below. If b k dominates over the remaining bj-s. j = 
0 ,,k — 1 , the integration of the dynamical system 


(b k ) s + ( 2 k + j^[) bb k + 2( |“ og ^| fc - 0 
b s + 2 6(0 — b) = 0 
a = b k (i + i^r) 

ds_ _ 1 _ As _ „ 


( 1 . 20 ) 


leads to finite time melting with the rate (II.4p for k = 0 and (11.51) for k > 1, which 
as solutions to (ll. 20 p display k instability directions. 


step 3 Energy estimate. We now construct a solution of the form 

k 

v(s,y) = ^2bj(s)r] b ^j(y) + s(s,y) 


j =0 


with 

( £ ,Vb,j)=0, 1 <j<k (1.21) 

and close an energy estimate for the remainder s. Here we note that for a solution 
to the linear problem 

d s e + T-L b £ + (a — b) Ae = 0, 

the time dependance of b(s) yields a modified energy identity 

J £ 2 e dy = ~(H b £,£) + (b s + 2b(a - b))^ J |y| 2 e 2 e '^ L dy 

and hence the choice of b (11.191) to cancel the second term in the energy identitjjj]. 
To the leading order, thanks to the orthogonality conditions (11.211) and the spectral 


Avhich involves a different type of norm for which the spectral gap constant is not explicit and 
would thus lead to severe difficulties. 
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gap estimate in weighted spaces associated to the knowledge of the kernel of Tib, we 
obtain the fundamental energy estimate: 

Id f o b \y \ 2 f o b \y \ 2 

J e 2 e 2 dy = -(H b e,e) < -{2k + 2)6 j e 2 e 2 dy. 

An integration-in-time will produce the necessary decay to close the bound on the 
radiative part of the solution, i.e. e. The situation is however more complicated 
since the problem cannot close at the level of H 1 Sobolev regularity, and instead 
forces us to take one more derivative. However, at the H 2 level the corresponding 
energy identity produces dangerous boundary terms. These come with a particular 
structure and may beautifully enough be handled through time integratior@, see 
Proposition 13.101 This part of the analysis is a replacement for the polynomially 
weighted estimates in [42] . and uses in an optimal way the dissipative structure of 
the equatior@ and the nonlinear algebra induced by the free boundary. 

The construction of the manifold of initial data to avoid the codimension k instabili¬ 
ties of the system of ODE’s (11.181 is finally performed using a now classical Brouwer 
type argument as in mmmm- 

step 4 Freezing. These regimes correspond to an expansion of the circular ice block, 
reflected in a change of sign in (11.161 . (11.171 : 

^ = A>0. 

A 

This causes a modification in the spectrum of the operator 
Hb = —A — BA, v{\fB) = 0, B > 0, 
which admits shifted eigenvalues: 

A S ,*= 24 + 2+1^ + »(|j^j). 

Computations parallel to the melting case lead to the dynamical system 

(B k ) s + (2k + BB k + 2( f log A if fc = 0 

B s + 2 B(B - A) = 0 
A = B k (l + 

ds _ 1 \s_ _ A 

K dt ~ \ ~ 

which after time integration produces the global-in-time freezing regimes (11.81 . The 
energy method is run along similar lines for very well localized initial data, since it 

By 2 

now involves the confining measure e 2 dy, B > 0. The analysis is slightly simpler 
thanks to the better decay of the B k mode which induces a stronger decoupling 
from the remaining modes. 

1.5. Plan of the paper. In section [2] we use a Lyapounov-Schmidt like argument 
to compute the bound state of and the associated spectral gap in weighted norms 
in both the melting regime 4*- < 0, Lemma 12.61 and the freezing regime 4^- > 0, 
Lemma 12.81 In section [3] we construct the melting regimes. We introduce the 
nonlinear decomposition of the flow, section 13.11 compute the modulation equations 
using the free boundary geometry, sections 13.21 and 13.31 and close the key energy 

8 This is reminiscent of similar essential issues in [32 :. 

^whereas the energy method in □IllIQdD works in both the dissipative and dispersive settings, 
but barely uses the dissipative terms in the energy estimates. 
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bound, Proposition 13.21 The proof of Theorem 11.11 now follows from a classical 
shooting argument a la Brouwer detailed in section 13.61 In section [4j we deliberately 
follow a parallel plan for the construction of the global-in-time freezing regimes. 


2. Spectral analysis in weighted spaces 

We compute in this section the k first eigenvalues of the linear operator 

?4,± = — A bA with boundary condition u(l) =0 

and the associated spectral gap estimate in the perturbative regime 0 < b < b*(k). 
The proof relies on a Lyapounov Schmidt type bifurcation argument at b = 0 per¬ 
formed in weighted Sobolev spaces. 


To ease notations, we fix 

PI 2 

± = —, 'Hb = 'Hb-, p = p~ = e 2 , b> 0, in sections 12.21 f2~3l 12.41 (2.1) 

and we omit the — subscript for the sake of simplicity. The case b < 0 with the p + 
weight and the operator T~Lb,+ is addressed in section 1231 

2.1. Coercivity for the harmonic oscillator. We recall in this section without 
proof the classical estimates for the harmonic oscillator. 


Melting case: consider —A + A on (H^ _ 0 , (•, -)o). This operator is self adjoint for 
the (•, -)o,- scalar product as is easily seen by writing 

—A + A =-— d z ( p-zd z ). (2.2) 

P-Z 

The normalized Laguerre polynomials |45j 

p x pk 

L k (x) = -^(e- x X k) , keN, (2.3) 

solve 

XL'l + (1 - X)L' k + kL k = 0 

and hence 

P fc (r) = L fc (y) (2.4) 

diagonalizes the harmonic oscillator: 


(—A + A )P k = 2 kP k , (P n , P m ) p ,o = 1- 
Moreover, they satisfy the double normalization condition: 


r+oo 

/ L/ n L m e dx — d nm , L n (t)) — 1 
Jo 


or equivalently 

(PjiPk) o,- = Sjk, Pk{ 0) = 1, 

and the classical induction formula 

A P k = 2 k(P k - P k _ i), k > 1. 


(2.5) 


( 2 . 6 ) 

(2.7) 


An extensive overview of Laguerre polynomials can be found in |45j . We recall the 
standard sharp Poincare inequality for the harmonic oscillator: Vit G Hj } _ 0 with 

(u,Pj) o,-= 0, 0 <j<k, 
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there holds: 

\\9 z u\\l 2 > (2k + 2 )|| iz ||^2 ■ 

P— >0 P-,0 

Freezing case: Consider —A — A on (Hj 0 , (•, -)o). Then the map 


( 2 . 8 ) 


T 2 


T 2 
0 ’-2 


u i —y w — e 2 u 

is an isometry and integrating by parts: 


[ 9 M 2 

/ |Vu| e 2 p + dz = 


/ \Vw\ 2 p-dz + 2B / 


or equivalently: 

w 2 

(—A — A)v = (—Aw + Aw + 2w) e 2 
Hence the family of eigenvectors 


(2.9) 

( 2 . 10 ) 


i*r 

Pj = Pj e ~ 


diagonalizes the operator, and there holds the spectral gap estimate: Vrt £ Hj J+ 0 
with 

(u, Pj) o,+ = 0, 0 <j<k, 

there holds: 


'M\l 2 > (2A: + 4)||u||^ 2 

P+,0 ^ O 


P+,0 


( 2 . 11 ) 


Remark 2.1. The shift 2 = d in (12.1111 . where d stands for the dimension of the 
ambient space, will be crucial for the computation of the freezing rates. 

2.2. Almost invertibility of the renormalized operator. Recall the notational 
convention (ED. We consider the renormalized operator 

Hi, = —A + A with boundary condition u(Vb) = 0 

in the radial sector and for 0 < b < b* small enough. Thanks to the boundary 
condition u(sfb) = 0 and (12.2|) . H\, is self adjoint for the scalar product (■,•)& on 
the Hilbert space Hj jb given by (11.111) . We claim a near invertibility property of Hi, 
which is the starting point of the Lyapounov Schmidt argument. 

Before stating the lemma, we introduce some notations. We first fix a frequency 
size 

A'eN 

and a sufficiently small parameter 

0 < b < b*(K) < 1. 

Universal constants in the sequel may depend on I\, but are uniform in b € 
(0,6* (K)). We define the Gramm matrix 

M b,k = ({Pi,Pj)b)o<i,j<k, k < K. (2.12) 

Observe from (12.51) that 


(Pi, Pj)b = (Pi, Pj)o + 0(1 zdz) = 6 U + 0(6) 

J\z\<Vb 


\z\<Vb 

and hence 

Mb.k = Id, + 0(6) is invertible 
for 0 < 6 < b*(k) small enough. We introduce the vector: 

Pk (Pj)0<j<k 


(2.13) 

(2.14) 
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and consider the function 

m k (b,z) = (M^Vk(Vb),V k (z)), 

which by (12.311 and (I2.14|) satisfies: 

k 

m k (b, z) = J2l 1 + °( b )] p j( z )- 

j =o 

We now claim: 


(2.15) 


(2.16) 


(2.18) 


Lemma 2.2 (Near inversion of H 5 — 2k). Let k € N and 0 < b < b*(k) small 
enough. Then for all f € L 2 b with 

(f,Pj) b = 0 , 0 <j<k, (2.17) 

there is a unique solution u € Hi, to: 

Pi° 

H b ,kU = / where Hb,kU = (Hb — 2 k)u — Vbmk{b , z)d z u(Vb) 

(u,Pj) b = 0 , 0 <j<k. 

Moreover, 

||Au || i2 + \\d z u \\ L 2 + ||Au || L 2 + ||u || i2 + |log&||\/&0 2 u(\/&)| < ||/]| L 2 (2.19) 

p,o p,o p,o p,o 1 1 p,o 

Proof of Lemma \2.2\ We use a Lax Milgram type argument in H p b . Let k € N and 
define the constraint set 

C := {u€ H X ft b | (u, Pj) h = 0, 0 < j < A:}. 

We consider the problem of minimizing the functional P : H p b —>• M. : 

P{u) = / |cLu| 2 /ozcfo — 2 fc / u 2 pzdz — (f,u)b 
J z>s/b J z>\/b 

over the constraint set u € C. Let 

A, = inf ^(u). 
uSC 

We recall from the standard Poincare inequality for the harmonic oscillator ()2.8[1 
and the compactness estimate (IA.1I) that for spherically symmetric v € H p 0 with 
(v , Pj)o = 0, 0 < j < k, there holds: 

j \d z v\ 2 pzdz — 2k J v 2 pzdz> J (1 + \z\ 2 )\v\ 2 pzdz. ( 2 . 20 ) 

Applying this to v(z) = ul„ >v ^ € H p 0 , we conclude that lb > —00 and that 
any minimizing sequence u n is uniformly bounded in H pb . Therefore, up to a subse¬ 
quence, using the compact Sobolev embedding H p b <-^ L 2 b and (12. 2011 . we conclude: 

u n ^ u in H p b , u n ^ u in L 2 b . 

In particular, using the local compactness of the embedding 77 1 (M) C L°°(R), this 
implies 

u(Vb) = 0 , (u, Pj)b = 0 , 0 < j < k 

and u is a minimizer of T over C. By a standard variational argument, there exist 
Lagrange multipliers A j € M such that 

k 

H b u = f~J2 X > P r 

3 =0 


( 2 . 21 ) 
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Hence from standard regularity argument, u € H 2 oc (r > y/b). We may then take 
the scalar product with Pi and compute: 

(H b u,Pi) b = -f —d z (zpd z u)Pipzdz = Vbe~ b/2 Pj(Vb)d z u(Vb) + {u,H b Pj) b 
Jz>Vb Zp 

= Vbe~ b / 2 Pi(Vb)d z u(Vb) + 2 i(u,Pi) b = Vbe~ h ^ 2 Pi(Vb)d z u(Vb). 


We conclude from (12.211) . (I2.17p . (|2.12|) that 


Vbe b/2 d z u{Vb ) (Pi(Vb)) = -M 6 >A .(Aj)o<j<fc 

V / 0<i<k 

or equivalently 

(Aj)o <i<k = -Vbe h/2 d z u{Vb)M b lVk{Vb) 

and hence u solves (12.181) from the definition (12.151) and (12.211) . We now observe 
that u € C ensures 

(■ m k (b,-),u) b = 0 

and hence taking the scalar product of (|2.18l) with u, using ( 12 . 201 ) with v = ul 2>v ^, 


and the identity (H b u,u) b = \\d z u \\ 2 r2 yields: 

L P ,b 


W u Wl\ £ \\dzu \\ 2 L 2 - 2 k\\u \\ 2 L 2 = ( f,u) b 

p,b p,b p,b 


and hence 

IMI l\ + \\dzu\\ L 2 < ||/|| L 2 

p,b p,b p,b 

We now integrate by parts to compute: 

( 2 . 22 ) 

(H b u,logz) b = (u, l) b - ^\logb\Vbe~ b/2 d z u(Vb). 

(2.23) 

We estimate from (12.161) 

ll"dc(Mllz, 2 , Z 1 

p,b 

and hence (12.281). (12.181) ensure 

(2.24) 

\^d,u(Vb)\ < |log(>| (ll^Hli,;,, + IMI 14 ,) S , logt | (ll/llil, + 

| Vbd z u(Vb)\j 


which together with (I2.22p yields 


\\d z u\\ L 2 b + |M| L 2 pb + \logb\\Vbd z u(Vb)\ < 

We now use the equation (12.18P again and teh bound (12.241) which yield 


L 2 ,- 

P,b 


\\H b u\\ L2 <\\f\\ L2 

p,b 


We then use a Pohozaev type integration by parts to compute: 

f zd z u(H b u)zpdz = — f d z (zpd z u )zpd z u — 

J z>y/b J z>y/b P 


z>Vb 

1 1 1 +0 ° 

-Z p (d z u) - 
2 PiVb 


\ [ z 2 p 2 {d z u) 2 ^V-dz 

2 Jz>Vb P 

he~ b l 2 \d z u{Vb )\ 2 + \ [ (A ufzpdz 

2 2 Jz>Vb 


= o 


2 

Ll , 


+ 




zpdz. 


which after a simple application of the Cauchy-Schwarz inequality yields ||Au || L 2 


r 2 and (12.191) is proven. 

p,b 


p,b 


□ IA 
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2.3. Partial diagonalization of H b . We are now in position to diagonalize Hi, for 
frequencies 0 < k < K under the smallness 0 < b < b*{K). 

Proposition 2.3 (Eigenvalues for H b ). Let K € N. Then for all 0 < b < b*(I\) 
small enough, H b admits a sequence of eigenvalues 


Hbifb,k = Xb, k i/j b ,k, ifb,k € H ] pb , 0 < k < K, 

such that for each 0 < k < K, the following properties hold: 

(i) Expansion of eigenvalues: there holds the expansion of the eigenvalue 


A b,k 2/c + , . + A & A b k O 

|log&| 


|log6| 2 ) ’ \^ b ^ b ’ k \ ~ 6|l 0 g6| 2 ' 


(ii) Expansion of eigenvectors: there holds the expansion 
j ifb,k = T btk (z) + ifb,k{z) 

} T b ,k{z) = P k {z) logz-log(Vb) + YljZo hb,jkPj (z) logz-log(Vb) 


(2.25) 


(2.26) 


(2.27) 


with 
f^bjk = 


(k-j)\logb\ 


+Tb,jk, hb,jk — ° (| log £| 2 ) ’ - O (|i og6 | 3 ) ( 2 - 28 ) 


and 


||Al/: 6ifc || L 2 + ||Z lfb,k\\L 2 , + HV+fcll// 1 , + ||Al/ ; 6,fc||L 2 , + ||^feV ; 6,fc|lz, 2 , 

p,b p,b p,b p,b p,b 

+ | log&| | Vbd z ^b, k (Vb) | < — 

\\d b ifb,k\\H\ + \\ogb\\Vbd b d z ifb,k{Vb)\ < 


and 


IIA^felL 2 , + \\z 2 ifb,k\\ L 2 h ~ l lo £ 

p,b p,b 

(iii) Spectral gap estimate: let u € H pb with 


(2.29) 

(2.30) 

(2.31) 


{u,ifb,j)b = 0, 0 <j<k, 


then 


Moreover, 


%u\\ L \> 

p,b 


2k+ 2 + 0 


|logb| 


\ u W L 2 ■ 

p,b 


\b,o = inf 


p,b 


6^,6 Mlh 


p,b 


and 


ifb,o(z) > 0 for z > Vb. 

(iv) Further identities: there hold the algebraic identities for 0 < j < k: 


|logb| 


2 r 


(' lfb,ki lfb,k)b 

{bd+fb.j i lfb,k)b 1 

(ifb,k,ifb,k)b |logb| 


1 + 0 


|logb| 


fijk + O 


|log6| s 


2 / ’ 


(2.32) 

(2.33) 

(2.34) 

(2.35) 

(2.36) 
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Moreover, 


IIA ^b,k - 2fc(' lfb,k - '06,fe-l)||lf2 < 1, 

p,b 

1,1 ^ + 5 lip 

ll" 2 ^,^- + (2/c + 2)'05 j fc_)_i — (4/c + 2)'06,fc + 2ki/jb,k-i \\h 2 p b 1$ 1- 


(2.37) 

(2.38) 

(2.39) 


Remark 2.4. From standard Sturm Liouville oscillation argument, if b . k vanishes k 
times onz> y/b, and hence only the ground state if b) o is nonnegative. 

Remark 2.5. Constants in Lemma f2. 31 depend on the frequency K but are uniform 
in 0 < b < b*(K). 


Proof of Proposition \2.iA The argument is of Lyapounov-Schmidt type. We remove 
the b subscript as much as possible to simplify notations. 


step 1 The Lyapunov-Schmidt argument. Let T k {z) € Hj jb be given by (|2.27p . and 
introduce the universal profile 

Tb(z) : = log 2 - log (Vb). 

Then 

J (Hi - 2 k)T k (z) = Q t (z) + Sjto 1 Ht [20 - k)T b (z)P j (z) + Q,} 

l Qj = -t^f- + Pj(z)- 

Observe from (12.4[i that Qj = Qj(r 2 ) with degQj < j — 1 and hence 

VO <j< k, Qj € Span(R ? )o<j<fc. 


(2.40) 


(2.41) 

(2.42) 

(2.43) 


We solve the eigenvalue problem 

(■ H b - 2k)if = n k if, 

by representing if in the form 

if = T k +if. 

and hence (|2.42jh (12.43|i give an equation for if : 

(H b -2k)if = -(H b -2k)T k +/j, k (T k + if) 

k -1 

= -Qk(z) ~ Ajfc t 2 (i - k)%(z)Pj(z) + Qj] + fi k (T k + if) 
j =o 

= : f(if). (2.44) 

We define {fij(if))o<j<k by imposing the relations: 

^p- h J} )b = Vbd z if(Vb ) (M~lV k {Vb)) , i = 0,...,k, (2.45) 

which is proved below to correspond to an invertible linear system on ((AL',fc)o<j<fc—l> /•*&)■ 
Observe that (|2.45p allows us to rewrite (12.441) as: 

k If PA. 

(2.46) 


=i? w 


3=0 "■ 


with H bk given by (|2.18p . Thus to find if, by Lemma 12.21 we are left with solving 
the fixed point equation 
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and we indeed claim that the operator H b } o F is a strict contraction from the 
closed ball 


B n := 


€ H} tb 


IIA-0II L 2 +Vb\d z i>(Vb)\ + 11-011 H\ < ppyr 

P’ b p* |logo| 

ip(Vb) = 0, ('4>,P j ) b = 0, 0 < j < fe|, 


(2.47) 


to itself for a universal large enough. 

Computation of fijk, //&: We invert (|2.45|> . Indeed, we rewrite 


/( 0 ) 


-Qfc + 


i + 


21ogz 

|log6| 


+ /bcV’ 


E 


3=0 


|log6| 

2 


[2(j - fc) - «fe] Pj 


1 + 


21 og 2 

|log6| 


+ 


|logb| / ' 


Let 

M (A^OAu • • • > f^(k— l)fc> A^fc) an d |/l|oo • l/U'fcl 0<j<k —1 j | t^k I} • 

Using the almost orthogonality (12.131) and the (■, -^-orthogonality of if and Pj, 
j = 0,... k, we obtain: 


{m,Pj)b 

(Pj’Pj)b 


{f(4>),Pk)b 
(Pk, Pk)b 


k -1 


(fc - j)|log6|//jfc - {Qk,Pj) o + (0/2) j + /i fc y^ y fi£k{Pj,PeTb)b + 0(|b|), 

r=o 

(2.48) 

Ao J ^ - (Qk,Pk) o + Mfc ^ Hek(Pk, PfPb)b + (Cfik + 0(|6|), 

£=0 

(2.49) 


where O = (U\ ? )ij=o....A: = 0(1), is a (/c + 1) x (fc + l)-matrix with bounded entries 
in the regime where b is small and (0/2)j = Oj£/t£fc + CikUk is the 2-th entry 

of the vector 0/2, i = 0,..., k. The above system contains quadratic terms and 
it can be solved for /2 by the following simple iteration argument. Given ft = 
(Aofc, ■ ■ • -/2(fe-i)fe, Wfe)> consider the system 

= (* - j)|log6|/i ilfc - <Q fc , Pj)o + (C'/2)j + 0(|6|), 

\Pj i Pj/b 

{ -§^ - + «?M) t + 0(\b\), 

\-tfe) *k)b Z 

where Qj = O/j + XJEq fak{Pj, PeT b ) b , j = 0,... k- 1. Assuming that |/2|oo < 1, 
we have O = 0(1) and we can invert the above system to obtain 


ftjk 


(Qk,Pj)o ( 1 y/b\dMVb)\ \ 

(k-j)\logb\ y |log6| 2 |log6| J 


,0 <j<k 


1, (2.50) 


2(Qk,Pk)o ( 1 Vb\dMVb)\ \ 

|log6| y |logfc>| 2 |log6| J 


(2.51) 


and therefore \fI\oo ]k^b[ — 1 f° r b < b* sufficienty small, where we used (12.471) . 
Contractive property follows easily in a similar manner and for a given ip € B a we 
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obtain the unique solution jl. We may moreover integrate by parts and use (12.61) to 
compute: 

(Qk,Pk) o = (—2^ + P k ,Pk) o = P*(0) = 1. 

Similarity for j < k — 1: 

(Qk,Pj) o = -2{^,Pj) o = 2 + (^,P k )o = 2 

p> 

since -f is a polynomial of r 2 of degree < j — 1 < k. Note that by (12.501) and (12.511) 
we additionally have the bound 


\u\ 


OO 


|log6| 


+ 


\Zb\d z 4>(Vb)\ 

|log6| 


(2.52) 


Estimating Let us introduce the approximate projection operator: 


IP,/ := E 


(f,Pj)b „ 

'P; ll 2 r 


1=0 


Observe that (I — P k)Pj = 0(6) and by (I2.41|) (I — P k)Qj = 0(6) for any j € 
{0, 1,... , k} and the almost orthogonality relation (12.131) . From (12.471) we have that 
(I — Pfc)^ = 4>. Thus, from (12.441) and (12.46}) we obtain that 

|log6| 


F W = - IP,) 


P fc logz + P k 


+ Pki’ 


k -1 


+ Pj,k [2(fe — j) + Mfc] (I — Pfc) 


l=o 


Pjlogz + Pj 


|log6| 


+ 0(6(1 + |/f|oo)) 


k -1 


: Mfc(I - Pfc) [-Pfclogz] + + E i 2 ( k ~ j) + Atfc] (! - F fc) [i^Tog^] 

1=0 

+ O (6(1 + l/Zloo) + 6|log6||/x| 00 ). 


Therefore 


\\nn\L\ 

p,b 


\iAc 


i + 


IT 2 
' L P,b 


< c 


1 + \Z6|5 2 ^(\/6)| 

|log6| 


+ 0(6(1 + |/7| c 

1 + 


+ 6|log6p| c 


a 


|log6| 


(2.53) 


for 0 < 6 < b* with 6* sufficiently small and a universal constant C > 0. Applying 
Lemma [2.21 and using (12.531) we conclude from (I2.46|) 


II a ^IIl 2 + Ill’ll h\ + |log6||v / 63 z i/((v / 6)| < ||F (^)|| L 2 

p,b p,b p,b 

1 + Vb\d z ii>(Vb)\ f a \ 

|log6| V |log6| y 

and hence 

ll A ^l \l 2 ,+ IIV’II h\ + |log6||\/6<9^(\ / 6)l + \\PWWl 2 , < pAr 

P,b P,b p,b |iOgO| 


and ijj € B a for a > 0 universal large enough and 0 < 6 < b*(k) small enough. 
Therefore ij) € B a . 
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Contractive property. To show the contractive property, note that for any 0i, 02 € 
B a by (12.4411 we have that 


/0 i) - /0 2) 
|log&| 


2 

k -1 

E 2 

i=o 


(/bcOi) ~ Pk(h)) p k 
|log6| 


1 + 


21ogz 

|log6| 


+ (/hfc (0 1) - Wfc( 0 2)) 01 + Wfc(02) (01 - 02) 


(Mjfc(0 1 ) - fJ-jkih)) 1 [2(j - fc) - // fc ] Pj 


1 + 


21ogz 

|log6| 


+ 


2Q y 

|log6| 


By a calculation analogous to (12.4811 and ()2.491) . we can evaluate the (•,-^-inner 
product of f (0i) — /(0 2 ) with Pj, j = 0,1,..., k and thereby estimate |/n/fc(0i) — 
/ijfc(0 2 )|, |/ifc(0i) - Mfc(02)|- Using (|2.47|) . we arrive at 

ll/(0l) - /(02)||x,3 < 7r~TT ||01 - 021|£, 2 , 

P . 6 | logo) P.& 

which together with F = (ld-F k )f gives ||F(0i) - F(0 2 )|| L 2 & fc\\ L 2 pb - 

The operator H b l is continuous by Lemma [2.21 and therefore for a sufficiently small 
0 < b < 6* the operator H^ 1 o F is a strict contraction. By the Banach fixed point 
theorem, there exists a unique 0^ € satisfying (12.251) . The Frechet differen¬ 
tiability of 0b with respect to b at any fixed b > 0 follows similarly, the classical 
details are left to the reader. 


step 2 Proof of (12.311) . We estimate from (I2.28p . (12.271) : 

l|ATb, fc || i2 < |log6| 

p,b 

and (12.291) now implies 


||A0& fc|| L 2 , 

p,b 


;$ |iog&|. 


We may now apply (IA.1D to zipb and (12.311) follows from (12.541) . 

step 3 Spectral gap estimate. Let u € H * b satisfy 

(u,0bj)b = 0, 0 <j<k. 

Let 


v = 111 


z>y/b 


E 


('ll, Pj)b 


\Pj r 2 

j = 0 11 ■> '' h p ,0 


Pi € H 


P, 0 


then by the Poincare inequality 


ll^wll l 2 > (2A: + 2)||u|| L 2 . 

p,b p,b 

We now compute from (12.271) for 0<j<k: 

/ j ~ l ' 
0 = (u, ipj) b = ( u, Pj(z)Tb{z) + y Pi(z)-Pi(Vb ) T 06,fc 


i=0 


and hence using (12.281) . ([2.29)1 : 

l(“> p j)b| 0) 


1 


|log6| 


«Hl2 0 <j<k 


(2.54) 


(2.55) 


p,b 


< 


\u\ 


L 


2 

p,b 


|log6| 


from which 
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Injecting this into (12.551) implies 


M\ L 2 > 

p,b 


2k+ 2 + 0 


|log6| 


p,b 


(2.56) 


To prove (12.331) . (|2.34D . let 


Pb = inf 
u£C b 


\\9 z u\\ 2 L 2 


p,b 


\ u \\ 2 t2 


p,b 


Then any minimizing sequence normalized by ||ttnHz/ 2 = 1 is bounded in H^ oc (r > 

Vb). By the compactness of radial Sobolev and trace embeddings and (IA.1I) it 
strongly converges in L 2 b . Hence the infimum is attained and from Lagrange mul¬ 
tiplier theory: 

Hb4>b = Pb4>b- 

Moreover, since \(pb\ is also an infimum, we may assume (j) b > 0. If p b 7^ A&, then 
(<^6)'06,o)6 = 0 and hence from (12.321) : 


Pb = WdzMl^ > 1- 

p,b 

\\9z^b\\ 2 L 2 

Note that A b,o = ppp p ' b > Pb by the definition of pb . Together with (12.261) this 

L p,b 

contradicts the definition of pb for 0 < b < b* small enough. Hence A^o = Pb is 
the bound state. The simplicity of the first eigenvalue follows from a classical ar¬ 
gument [6]- and hence ^ = 4>b > 0. Note that ^ > 0 for z > \fb by the strong 
maximum principle. 


step 4 Estimate for db\,k- Applying d b to (I2.25p . we obtain 

Hbdb'l/jb,k dbXbjk^bjk T A b,k^b^b,k- (2.57) 

Evaluating the (•, •)& inner product of (12.571) with anc l integrating by parts we 
obtain 


^bXb,k ||^6,fe || t 2 T Xb } k{db^b,ki i>b,k)b — (Hb&b' , Pb,ki'^b,k)b 

p,b 

= {db^b,kiH b ^b,k)b ~ d h ^b+{Vb)d z ^ b +{Vb)p{Vb)\/b 
= A b ,k(db^b,k,A,k)b ~ d b ^b+{^)dzi> b +{Vb)p{Vb)Vb. 


Therefore 


d b X b,k = ~ 


d b MVb)d z MVb)p(Vb)Vb 


II 2 

Ml 2 . 

P,l 


(2.58) 


From ipb,k(Vb) = 0 it follows that <9 b^b,k{Vb) = ~-^^d z ipb,k(Vb) and therefore 
from (I2.58P 

\d z i’ b ,k(Vb)\ 2 p(Vb) 


&h^h.k 


2||'0t.,fc ||?2 

L p,b 


In particular, since |<9 z V’&,fc( 1 /&)| = 0(-^=) by (12.271) and (I2.29p . and 
|log&| 2 by (I2.28P - (I2.29P it follows that 


(2.59) 


112 > 

,/c 11 T 2 rsj 

P,b 


| d b X bt k | < 


1 

6|logt»| 2 ’ 


which is the last claim of (|2.26p . 


(2.60) 
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step 5 Estimate for \d b /j, b i k \, i = 0,. .., k — 1. From (12.271) we obtain 

^ k —1 ^ fc—1 

Pfc + ^ d b fi b j k PjTb - ^ y] Hb,jkPj- (2.61) 

l=o l=o 

From (I2.27p and (12.571) it follows that 

Hbd b ^b,k = Fk + \,kdbipb,k, (2.62) 

where 

-Pfc = -H b d b T b)k + d b (\ b)k T b)k ) + d b \ b} k^b,k- (2.63) 

Rewriting (12.621) in the form H b d b ^ b ^ = (d b X btk ^b,k + X b)k d b ^ b)k ) + A b) kd b T b) k and 
evaluating the (•, -^-inner product with Pj, j = 0,..., k — 1, we obtain 

(H b d b 'ip b:kl Pj) b = A bt k(d b T btk , Pj) b (2.64) 


since (ipb,ki Pj) b = (d b ij) b ^ kl Pj) b = 0. On the other hand, from (12.611) we have that 


( d b T b>k , Pj) b = XIj k + '^2 dbHb,ik i M j } 2 ^ Pj)b) ~ ^ ^2 ^b,ikM- 

o—r\ \ 

in i 4- ni 

'6|log6| 


_ 

i =0 A 

|log6| ^ ~ . 1 


fe-l 


j=0 


2 1 26 

^bl^b,ik^ji + (2.65) 


i=0 


where (cjj)ij=o,...,fc-i = 0(1) an< l we used Mji = 5ji + 0(6) for j = 0,..., k — 1 and 
the first two claims of (12.281) which have already been proven above. On the other 
hand, observe that by integration-by-parts and the orthogonality (<9f,V’fe,fc) Pj) b = 0 
we have 


(■ H b d b T b)k , Pj) b 


dM^Vb^jiVb) - d b T btk (Vb)d z Pj{Vb2) . (2.66) 


By (12.611) \d b T b)k (Vb)\ < 1/6 and 

( k —i \ k— 1 

d z P k {Vb) + ^2 Hb,jkd z Pj(Vb) ) + d b fi b , jk Pj(Vb), 

l=o / l=o 

which together with (12.661) leads to 

k-l 

(H b d b T b)k , = p(Vb) ^2 dbVb,jkPj(Vb) + 0(1), (2.67) 

l=o 

where we note that d z P k (^/b) = \fbL ' k (|) = 0(y/b) implying maxj = o r .. i / c \d z P k (-\/b)\ 

0(1). To see that L ' k (|) = 0(1) observe that by the definition (12.31) of the k-th La- 
guerre polynomial, it follows that L' k is a polynomial of degree k— 1. From (12.641) . (|2.65p . 
and (12.671) we conclude that 

] ^Y^d b pb,ik^ j i + c *jid b p b ,ik = 0( ft|iogfti ), c l* = 0(1), i, j = 0,. .. ,fc - 1. (2.68) 

The system (12.681) is invertible for 0 < 6 < 6* sufficiently small, and as a consequence 


sup \dbfj, bj ik\ < 
i=0,...,k—l 


1 

6|log6| 2 ’ 


(2.69) 


which completes the proof of (I2.28p . 
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step 5 Estimate for Recalling that (i/j b ,k, Pj)b = 0 , j = 0,1,... k, by 

the construction of ipb,k, we conclude that (dbipb,ki Pj)b = 0 since ^ b ,k{Vb) = 0. 
Moreover, the spectral gap estimate (12.811 with u = d b ^b,k^ -. > ^+d b ^b,k{Vb)l Q<z<yJ r b 
together with the bound A b k = 2,k + A b k imply 

2k + O(tt^t) 

\b,k\\ d bA,k\\l\ ^ 0 , — ~ — ° g i ^ \\dzd b ^ k \\ 2 L 2 + C\d b ^b,k{Vb)\ 2 . (2.70) 

p ’ b 2k + 2 + °\^\) p ’ b 


Evaluating the (•, -)b inner product of (|2.62p with d b ipb,kj integrating-by-parts, using 
Cauchy-Schwarz, and the Poincare-type inequality (12.701) we obtain 

\\d z d b A,k\\ 2 L 2 ~ ll-^fclli 2 + ^\d z d b 'tp b ,k(Vb)d b 'ip b:k (Vb)\ + \d b i’b,k{Vb)\ 2 . (2.71) 

p,b p,b 

From ’ipb,k(Vb) = 0 it follows after differentiating with respect to b that 

\d b M^)\ < ±=\d z ik tk (Vb)\ < (2.72) 

To estimate Vb\d z d b 'ipb,k(Vb)\ we note that 


0 = {dbi> b ,k, l)fc = (dbi’b,ki H b \ogz) b 
= (H b d b 'ipb,kAogz)b + e~ b/2 d b i>b,k{Vb) + ^\\ogb\Vbd z db$b,k{Vb)e~ b/2 . 
Therefore 


\VbdMb,k{.Vb)\ < + |A 6)fe |||^ 6)fe || L 2 

~ |i^6| + ]i^6f 



where we used (12.621) in the first estimate and the Poincare inequality (12.701) in 
the last. Using the Cauchy-Schwarz inequality and plugging this bound and (12.72|) 
into (12.711) we obtain 


\\dMb,k\\ L l b £ ll^fcll i? p<b + ( 2 - 73 ) 

From (12.611) and (I2.69p we conclude that II7), fell r 2 < j-. To estimate IIIE 2 note 

’ p,6 0 p,b 

that from (|2.40p we have the identity: 


k -1 ^ k-1 

H b d b T b) k = 2 kd b T b) k + ^ d b ^ b ,jk [2(j — k)T b Pj + Qj] + -r ^ Pb,jk{ k ~ j)P- 


i=o 


l=o 


and therefore using (I2.6ip and (12.261) we have 
-H h d b T bk + A b,kd b T b ,k = 0(-rj— -rr)d b T bjk 


|logb| 

k ~1 1 fc-1 

+ ^ ^ d b ^ b ,jk [2 (j — k)T b Pj + Qj] + — fk b ,jk(k — j)Pj 


l=o 


l=o 


Hence, using (12.691) and fl2.28|> : 


1 
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Thus, using the definition (12.63)) of F k , bounds (12.601) . (I2.29p . and the previous 
bound we obtain 

1 


I F, 


k\\L 2 


< 


+b ~ b|logb| 


+ \d b X b , k \ (|| T b ,fc|| L 2 pb + ll'0fo,fc|lz.2 b ) 


< - + 1 


b b|logb| 2 
Plugging this back into (12.731) we get 


|logb| + 


1 


< 


1 


|logb| 7 ~ log 


zd b tph,k\\L 2 


< 


1 


“p. 6 ~ 6|log6| 

and therefore, by the spectral gap estimate ()2.8|) . just like in (12.701) . we obtain 

1 


\\d b iJJk,b\\L 2 


< 


'p,b ~ 6|log6| 


and the proof of (12.291) is completed. 


step 6 Proof of (12.351) - (12.391) . 

Proof of (12.35)1 . We estimate from (12.271) , (12.291) : 

|logb| 2 


(V+fc) V+fc )b — 


(Pk,Pk)o + O 


|logb| 


(2.74) 


and (12.351) follows from the normalization 
Proof of (12.36)) . We compute from (12.611) , (12.281) , (|2.29)) : 


bd b if b)k + -Pk 


H 2 ~ |logb| 

p,b 


(2.75) 


and hence using (12.35)) : 

(hOhPb.j i f b .k)b _ 4 

(^b,k,i>b,k)b 


(k>g&) 2 

4 

(log&) 2 


1 + 0 


log& 


(- : y,V’6,i)b + 0( 1) 


--\\o g b\(P v Pk) b + 0(l) 


1 


+ o 


|logb| V |!og^l 


this is (I2.36p . 

Proof of (12.37P : We compute from (12.270 . (12.7P : 


k -1 


AV’fe.fc = X.P k T b + Pk + ^2 Pb,jkMPjPb) + 


3=0 


k -1 


= 2A:[Pfc - Pfc_i]76 + -Pfc + ^ Hb,jkMPj7b) + 

j=o 

2k(lp bt k V+fc— l) + 

where the remainder estimate 


II H 2 


< 1 


holds due to (12.281) . (12.291) and (I2.37p is proved. 
Proof of (12.38p . Note that 


fc-i 


Pk = 


|logb| 


ipb,k - i’b,k - Pklogz - ^ HbjkPjTb 


i =0 
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and therefore 



where we used (|2.75j) and (12.2811 . (12.2911 . This concludes the proof of (I2.38|) . 
Proof of (12.3911 . We have from (12.2711 . (12.2911 : 


Let = z 2 Pk, then 

(—A + A)<f>fc = (2 k + 2)<f>fc — 4 Pk — 4APfc = (2 k + 2)<f>fc — 4 P k - 8 k{Pk — Pk- 1) 

and hence the relation 

2k{$k,Pk)o = (2fc + 2)(<& fc , P k ) 0 — 4 - 8k ie ($ k , P k )o = 2 + 4/c 
(2/c — 2)(4>j, Pfc_i) 0 = (2k + 2)(4?^, Pfc_i)o + 8/c ie ($*;, Pfc)o = — 2/c 
($k,Pj) o = 0, 0 < j < /c — 2. 

Since P k is a polynomial we conclude that there exists a c k E R such that 

= CfcPfc+i + (4/c + 2)P k — 2kPk-i- 

Since P k ( 0) = 1, we obtain by plugging in z = 0 into the above relationship: 

z 2 Pk = —2(k + l)P k +\ + (4 k + 2 )P k — 2kPh—i, 
which yields (12.3911 . □ 

2.4. Diagonalization of Tib- We are now position to derive the bound state and 
the spectral gap estimate for the operator PLb- 

Lemma 2.6 (Renormalized eigenfunction). Let K E N. Then for all 0 < b < b*(K) 
small enough, the renormalized operator 

PLb = —A + bA with boundary value u(l) = 0 

has a family of eigenstates ry, k satisfying: 


PbVb,k = b\ b ,kVb,k, 0 <k < K, 


(2.76) 


with the following properties: 

(i) Structure of the eigenmodes: there holds the expansion 



(2.77) 


with 
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(ii) Further estimates on the eigenvector: there holds 


I \yVb,k\\b V 


|log6| 


J 2r lb,k\\b ~ 


[logfr| 

by/b 


II bd b rj b: 


fc life V 


|logb| 

Vb ' 


(2.79) 

(2.80) 


Moreover: 


IIA Vb,k - 2 k(rj b>k ~ rj b ,k-i)\\b + -^=\\d y [A - 2k(y byk - rj byk -i)\ life 
+^\\kL b [Arj byk - 2 k(y byk - 77fo,fc—i)] ||b < 


(2.81) 


and 


2 1 2 
||26<9b?76j - Ar/ 6)j + J^VbjWb + -^\\d y [2bd h r] byj ~ Arj byj + -^^Vb,j\\b 


1, 


H-^ll'Hfe 


2 bd b r) byj - Arj b)j + 
(iii) Normalization: 

(Vb,k, Vb,k)b 


|log6| 

|logfc| 2 

4 b 


hb,j 


< 


1 


Vb\logb\ 


1 + 0 


|log6|y_ 


(2.82) 


(2.83) 


Remark 2.7. Observe that (12.811) . (|2.82l) . (I2.83P imply the bound: 

\\bd b r) byk - k(r] byk - y byk -i)\\ b + —^\\d y [bd b r] byk - k(r) byk - r] byk -{j]\\ b 

Vb 

+^||%6 [bd b rj byk - k(r] byk - r] byk - 1 )] || b < -^=. (2.84) 


Proof of Lemma \2.(A Given u : 0 —>• M, let v(y ) = u(Vby), it is straightforward to 
check that 

n b v = b(H b u)(Vby), \\d^v\\ b = b~V \\d e z u \\ L 2 , t € N. (2.85) 

We therefore let 

flb,k{y) = Vb,k( z ), z = yVb. (2.86) 

and (12.761) follows. The estimate (12.811) follows by rescaling (12.371) . and (12.831) by 
rescaling (12.351) . The decomposition (12.771) follows from (I2.27J) . and (12.291) implies 
(12.781) . The bounds (12.791) follow by rescaling the bounds: 

\\ z Vb,k\\l? h + \\ z2 Vb,k\\L 2 h < l lo g fe l- 

p,b p,b 

Directly from the definition (I2.86|) . we compute: 


bd b r] byk 


Ay byk 

2 


+ bd b i/j b , k 


('/by), 


(2.87) 


which together with (I2.30p . (I2.3ip yields: 




|log6| 

V~b 


thus proving (12.801) . From (12.8711 : 


2 bd b rj byk - Arj byk = 2 [bd b Vb,k] (Vby) 
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and hence from (12.381) : 


2bd b r] b)k - Ag b)k + 


j ~Vb,k 


1 


1 


|logfc| 

and similarly for higher derivatives 


<^IW W +Hog 6 | 


i>b,k\\l£, ~ 


Pb ~ |logb|\/b 


□ 


2.5. Diagonalization of TLb- We now change sign and consider the operator for 
S > 0 

Bb = — A — SA with boundary value u(l) = 0, (2.88) 

which is a self adjoint operator on Hg + given by (11.121) . 

Lemma 2.8 (Renormalized eigenfunction). Let K £ N. Then for all 0 < B < 
B*(K) small enough, the renormalized operator 

TLb = —A — BA with boundary value 7/(1) = 0 

has a family of eigenstates 

B\y \ 2 ^ ~ ^ 

VB,k = e 2 T)B,k, TdBf)B,k = B\B,KhB,ki ^B,k = ^B,k + 2, 0 < k < K, (2.89) 

with r/B, k ,\B,k given by Lemma \2.(A Furthermore, there hold the following proper¬ 
ties: 

(i) Structure of the eigenmodes: there holds the expansion 


b lap 


VB,k = S B ,k(y)e 2 +t) 

SbAv) = p k{VBy)\ogy + Y!jZ] hB,jk p j(VBy)logy 


(2.90) 


with 


II A5 ?baIIl| 

Vb 

< —-— 

~ Ilog51’ 


+ \\dyfjB,k\\B + VB\\fj B ,k\\B + VB\\Afj B ,k\\B + \ \°ZB\\dyf)B,k(X)\ 

(2.91) 


(ii) Further estimates on the eigenvector: there holds 


\\yf)B,k\\b < 


| logS | 


\\Bd B f)B,k\\b 


B 

< | logS | 

~ Vb ' 


' 2 "Afcll B < 


\\y 2 VB,k 


[logs | 

bVb 


(2.92) 

(2.93) 


Moreover: 


+ 

+ 


\\BdBfjB,k - {k + 1 )\i)B,k+l — fjB,k]\\B 

=\\dy[Bd B fjB,k ~ (k + l)[f}B,k+\ ~ VB,k]\\B 

— ||Sb [Bd B fj B ,k - {k + l)[fjB,k+i - VB,k ]] lie ^ 


and 


2 1 2 
\\2Bd B fiB,j ~ A.fj B j + | Iogg | ??Bj||& + -^\\dy[2Bd B flB,j - A f) B ,j + j^"gj 




2Bdsf)B,j ~ ArjB,j + 


| logS | 


VB,j 


I B 


1 


\/S|logS| 


(2.94) 

VbjWb 

(2.95) 
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(iii) Normalization: 


( VB,k,VB,k)B = 


log B | 
4/i 


1 + 0 


log 51 


(2.96) 


Proof of Lemma \2.97[ This is a direct consequence of Lemma [2.61 Indeed, the map 


L b,+ -> L l,- 

b |„|2 is an isometry (2.97) 

v ^ w = e~2~ v 


and integrating by parts: 


/ |Vu| 2 e ^ p B ,+dy = / |Vw| 2 /9b _dy + 2B / |tc| 2 p s _dj/ 
or equivalently: 


(2.98) 


B\y\* 


PLbv = (—Aw + BAw + 2Bw) e 2 . 

Together with Lemma 12.61 this yields (12.891) . We now renormalize flA.ll) which 
yields: 

B t 5 l \\y k w\\B- < \\d y w\\ B - + Vb\\w\\ B -. (2.99) 

The isometric relation (I2.98|) and (12.991) imply the following comparison of norms: 

\\d y v\\ Bj+ < \\d y w\\ B - + Vb\\w\\b,+ (2.100) 

11 Ar; 11 = || A w - By 2 w\\ B - < ||Au;|| s _ + —7=\\d y w\\ B _ + 11^11^,(2.101) 


Vb' 


( 2 . 102 ) 


\\y k v\\ Bl+ < \\y k w\\ B -. 

Using (I2.99|) : 

I|Au||b,+ < ||Aw;||b+ B\\w\\b- + B\\Aw\\ b _ + B 2 \\y 2 w\\ B _ 

< ||Aw|| s _ + B\\w\\ B - + B\\Aw\\ B - + VB\\d y w\\ B -. (2.103) 
We also observe that 

\d y w(l)\ < \d y v(l)\ if w( 1) = 0. 

Using the above bounds together with (12.891) . (12.901) . (|2.78D . (I2.79p . and (12.831) yields 
(12.911) . (12.921) . and (|2.96l) . Moreover, 


\vf 


By* 


dBVB,k = \ d B r]B,k - 7f-VB,kj e 2 . 

Hence (12.931) follows from (12.801) . (12.791) . We now observe the fundamental conjuga¬ 
tion 

B\y \ 2 

h-VBj ~ 2 Bd B ffB,j = (Ar] B ,j — 2 Bd B riB,j) e 2 , 

and hence (|2.95|) follows from (|2.82l) . Moreover from (12.391) : 

| 2 „ 2 „ 


< _L 

~ Vb 


B\y\ rjB,k = z i’B,k = -(2k + 2)r] Bi k+i + (4:k + 2)ri B:k — 2kri B:k -i+F k , ||F fc || B 

and hence using (|2.84p : 

\v\ 2 1 

Bd B r] B ,k — B-^-r] B}k = k(ri B)k - rj B ,k-i) ~ ^ [ _ ( 2 ^ + 2)VB,k+i + ( 4 k + 2)rj B ^ k - 2kr) Btk _ 1 ] + F k 

= (k + 1) VlB,k+\ ~ 9 b, k] + F k with \\F k \\ B <-^= 


and similarity for higher derivatives, and (12.941) is proved. 


□ 
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3. Finite time melting regimes 


This section is devoted to the existence and stability of the melting process. In 
all the section, we let 

± = P = P~, b> 0. 

3.1. Renormalized equations and initialization. We start with the classical 
modulated nonlinear decomposition of the flow. We let 


u(t,x) = v(s,y), y = 


m 


A(0) = 1, 


(3.1) 


where A(0) = 1 is assumed without loss of generality thanks to the scaling symmetry 
(11.1411 . We define the renormalized time 

rt dr 

so > 1, 


s(t) = s 0 + [ 
■Jo 


A 2 (r) 


and obtain the renormalized flow: 


d s v + H a v = 0, a = -■ 


A ’ 


v(s, 1) = 0, d y v{s.. 1) = a. 

We now prepare our initial data in the following way: 

case k = 0 : We first claim that given 0 < b* <C 1 and e* with 

b* 


(3.2) 

(3.3) 


| e *|| 2 < 

1° 11 0 * rsj 


| log b* | 2 ’ 

there exists a locally unique decomposition 

b*r] b * fi + £* = br] bfi + e with (e,r) b ,o)b = 0, |6-6*|< 
Indeed, we define the map 

F(b, £*) = (b*r) h *,o - bri b;0 + £*,y b , 0 ) b 
which satisfies F(b*, 0) = 0 and 

1 log b * 12 r 


log b 


i* 12 ' 


(3.4) 


d b F(b*, 0) — — (rj b *,o + b*(d b r] b 0 )\ b=b * ,r] b * t0 ) b * — -- 


4 b* 


l + O 


log b* 


< 0 


from (|2.83p and the degeneracy (12.8411 for k = 0. The claim then follows from a 
standard application of the implicit function theorem. A Taylor expansion of F 
about (&,£*) = (b*,0) yields the bound 


I b-b* 

and hence 


log b : 


* 12 


46* 


l + O 


1 


log 6* | 
\b~b*\< 


< 1 
rs_/ I 

£ *\\b*\\Vb 

Vb* 

Vb* 


,o||b* ^ Ik* II b* 


log b* | 

2 Vb* 


which concludes the proof of (13.411 . 
We therefore pick an initial datum 


| log 6|* | log 6|* 


v o ~ b oVb* + £q, Ikollft* ~ | lo g 5*12 

and decompose the solution 

v(s,y) = 6 0 (s)% o(s ) >0 + e(s,y) with b(s) = b 0 (s) and (e, %, 0 ){, = 0 (3.5) 
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which makes sense as long as s(s,y) is small enough in the L 2 weighted sense. We 
let 


£2 — T~Li>£ 


(3.6) 


and define the energy 

£ ■= \\n b e\\l 


which is a coercive norm thanks to the orthogonality condition (13.5|) . see Appendix 
El We assume the initial smallness 


£(0) < b3{0) 

[ llogW 

(3.7) 

case k > 1: Let 



k + 1 
CM “ 2 k 2 

(k + 1 )a k 

, Ck, 2 = Ck,l ^ 

(3.8) 

Then we freeze the explicit value 



b(s) 

1 1 ^k,l 

2 ks slogs 

(3.9) 

and the sequence 

f K = 

l b 3 = 

k- f-1 | Ck,2 

2 ks ' slogs ’ 

o, 0 < j < k — 1, 

(3.10) 


where the index “e" stands for exact. 


Remark 3.1. Letting 


then (a e ,b, 6|) satisfy 


k + 1 


e . Cfc.i 

a = — -h 


2 ks slogs ’ 


(bl) s + bbl 2 k+ 2 


I log &l 


+ 


2(q e — 


I log b\ 


= O 


b, + 26(o' -6 ) = 0 

. - * (l + da) = 0 (^Lf) ■ 


s 2 (logs) 2 


We define 


k 

Qp(y) = ^2 b jVh,j(y) 

3=0 


and introduce the dynamical decomposition 


(3.11) 


(3.12) 


(3.13) 


v(s,y) = Qp( s ) + e{s,y), (r) b j{s), s) b = 0, 0 < j < k. 


(3.14) 


We let again 

£2 = 'HbS, £ ■= ll'HfeEll 2 , 

which due to the orthogonality conditions (13.141) is a coercive norm, see Appendix 
El We assume the initial smallness 


£(o)< 


fe 3 (0) 

|log6(0)|' 


(3.15) 


For k > 1, the set of initial data will be built as a codimension k manifold. To this 
end and in order to prepare the data, we consider the decomposition 


Ms) = bj(s) + bj(s), bj(s) = V ^ S \s , j = 0,... ,k. 

s(logs) 2 


(3.16) 
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Let the (2 x 2)-matrix A k be given by 

1 + 4 


j d k .— 


k(k + 1)' 


(3.17) 


The matrix A k is diagonalizable with one strictly positive )jl\ > 0 and one strictly 
negative eigenvalue //g < 0. Let P k be an orthogonal matrix diagonalizing A k , i.e. 


A k — A. k P k , A-k • — 

We define the new unknowns W k , W k -\ by setting 

■=Pk 


( Mi 0 

0 /i*. 


W k 

w fc _ 1 

and now assume the initial bound 

|W fe (0)| < 1 


| W fc _ 1 ( 0)| 2 + ^ 

j =0 


V k 

V k -! 


k 2 ^(0) 2 


< K l 


(3.18) 

(3.19) 

(3.20) 


for some universal constants K > 0, 0 < 5(k) < 1 to be chosen later. We then 
consider the bootstrap bounds 


£ < 


Db 3 

I log JI 


■j for k = 0, 
for k > 1 


- 1 M>L 

l |logb| 

for some large enough universal D = D{k) to be chosen later, and 
• for k = 0 

0 < bo(s) < b* 


(3.21) 


for k > 1, 
and 


|W*(a)| < K 


| W k - 1 ( s )\ 2 + ^2 

3=0 


k—2 T7- / \ 2 

v j( s ) 


< K 


(3.22) 

(3.23) 

(3.24) 


and define 


s = 


suo„^ { (13.2111 . (13.221 hold on [sq,s]} for k = 0, 


suPg> S0 { dSZH) ; (13.2311 . (14.13jl hold on [so,s]} for k > 1. 

The main ingredient of the proof of theorem 11.11 is the following: 

Proposition 3.2 (Bootstrap estimates on b and e). The following statements hold: 

1. Stable regime: for k = 0, s* = +oo. 

2. Unstable regime: for k > 1, there exist constants K, 5 = 5(K) <C 1 and 
(Vo(0), •••, 14 _ 2 ( 0 ), Wfc_i(0)) depending on e{ 0) satisfying (|3.15|> . (13.1911 and (13. 2011 . 
such that s* = +oo. 


Remark 3.3. Let us observe that our set of initial data is non empty and contains 
compactly supported arbitrarily small data in H 1 , see Appendix [G] 

Remark 3.4. The proof of the Proposition 13.21 is presented in section [331 
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From now on and for the rest of this section, we study the flow in the boostrap 
regime s € [so,s*). Note in particular the rough bounds 

\h\<b, \bj\ < —, 0 < j < A; — 1 and £ < — ^L= (3.25) 

I logo| Vllog^l 

for so > s 0 (K) large enough. 

3.2. Extraction of the leading order ODE’s driving the melting. We derive 
in this section the main dynamical constraint on the parameters (a, b, (^')o<i'<fe) 
which lead to the leading order modulation equations, and are a combination of the 
linear diagonalization of the PL b operator and the nonlinear boundary conditions. 


We start with the constraint induced by the boundary conditions. 
Lemma 3.5 (Boundary conditions). There holds: 


a = b i 

3=0 


1 + 


2a j 
|log6| 


+ o 


V£ \ 

|log6| 2 ' \logb\Vb)' 


+ 


e 2 (l) = ~a(a - b ), 


dyS 2 ( 1 ) — &s ^ 
3=0 


i 2a,- 

1 +- - 


|logb| 


+ 0 


b 2 


|log5|- 


(3.26) 

(3.27) 

(3.28) 


Proof of Lemma 1 3. ,71 We compute from (12.7711 , (12.781) , (12.2811 and recalling the def¬ 
inition (11.1311 : 


3 -1 


= 1 + Y,jf 


1 \ _ 2aj ( 1 

^ (j -*)|logb| ' ~ VMV + |log6| + V llog^l 2 


+ o 


i =0 


. (3.29) 


This implies that 

k k 

dyQhi 1 ) = bjdyrjbjj 1) = bj 

3=0 3=0 

Since v = Qp + e, it follows that 


1 + 


2a j 
|log6| 


+ o 


|log6|V ' 


£ y\ y =i ~ V V, y 


si dyQp\ v= i — \ dyQpi !) 


= a ~^Z b 3 

3=0 


2 aj 

1 +- - 


|log6| 


+ o 


|log£»| s 


2 / ’ 


which together with (1 A. 1011 yields (13.2611 . From (13.31) . v(s, 1) = 0 and d y v(s, 1) = 
— — = a‘ 

0 = TL a v( 1) = ( HbV + (a — 6)Au)(l) = e 2 (l) + a(a - 6), 
this is (|3.27p . Now from d y v{s , 1) = a, we have 

d s d y v(s, 1) = a s . 

On the other hand, taking d y of (13.311 . we have: 

0 = d s dyV + d y (PL b v + (a - b)Av) = d s d y v + d y s 2 + dyH-bQp + (a - b)yAv. 

We evaluate the above identity at y = 1. From (13.3p and d s v(s , 1) = 0, d y v(s, 1) = a: 

Aw(l) = aAu(l) = a 2 . 
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By construction, 


dyTibQp — 'y ^ ^b,jbbjdyVb,j i 

3=0 


and hence: 


+ d y £ 2 (i) + Y J Xbjbbj 

j=o 


i+ipri + a 2 {a — b) = O 
|log6| 


b 2 


|log6| 5 


We inject into the estimate the rough bound 

M ^ b 

which follows from (I3.25|) . (13.26|1 and (13.2811 is proved. 


(3.30) 

□ 


We now show how Qp is prepared to generate an approximate solution to (13.3|) 
with the suitable leading order dynamical system for (/3, A) induced by the spectral 
diagonalization of PL b - 

Proposition 3.6 (Leading order modulation equations). Under the a priori bounds 
of Proposition 1. 9. HI there holds 


d s Qp + PL a Q p — Mod + 'L 
where we defined the modulation vector 


Mod : = 


. ,, , 2(a — 6)6^ kbk ^ 

( bk)s + bbk\b,k H-n—7|-1— T ~^ 

|logb| b 


Vb,k 


(3.31) 

(3.32) 


fc-i 


+ E 

3=0 




VbJ j 


the deviation 

:= b s + 2 b(a — b ), 

and the remaining error satisfies the bound: 

ll *"‘ + 76 l|a “ 4 ’ llt+ s ||w ‘ 4,||t 5 + 

Proof of Proposition \3.(k By definition 

Pl a = Tdb + (a — b) A 

and we therefore compute from (13.1311 : 


(3.33) 

(3.34) 


dsQp( s ){y) T~^aQp — 'y ' 


3 = 0 


(bj)sVb,j + bs-^bd b r]b,j + bbjX b}j r] b j + (a - b)bjAr) bJ 


E 1 + bbj\b,j]rjb,j + (a - b)bj[Ar] btj - 2 bd b r] b fi + -j-bd b r]bp^ 

3=0 


E 

3 =0 


( bj) s + bbj X b j + 


2 ( a ~ b)bj ~\ jbfib _ 

|logb| \ Vb ’ J b [VbJ VbJ ~ l] 


+ (a - b)bj[Ar]b,j - 2 bd b 3j b j - rj b j] + ^<f> [bd b rj btj - j(r] b}j - r)b,j-i)\ 
The bounds (12.8211 . f|2.84|) . (I3.25|) . (13.301) now yield (13.3411 . 


□ 
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Remark 3.7. The presence of the |logb| in the denominator on the right-hand side 
of (13.341) makes it a true error term with respect to our bootstrap regime, and 
this term is one of the leading order errors when closing the energy estimates in 
sections 13.31 and I3~I1 


3.3. Modulation equations. From (13.31) . (13.311) we obtain the equation satisfied 
by the perturbation e : 

d s e + U a e = 7 (3.35) 

where 

7 = -Mod - T. (3.36) 

The nonlinear decomposition and the orthogonality conditions (|3.14|) generate a 
differential equation for the modulation vector (3 = (&j)o<j<fc in the setting of the 
bootstrap lemma [3721 which we now compute exactly. 


Lemma 3.8 (Modulation equations for bj ). 1. k = 0: the b law is given by: 

2 b 2 


b s + 


|logb| 


< 

~ |log6| 2 


(3.37) 


2. k > 1: the modulation dynamical system for the vector (6j)o<j<fc is given by 


(bk)s 4 — 
s 


bk + (k + 1) b] 
j =o 




~ 1 


+ 

( bk-i)s 4— 
s 


k-l 


J k -1 


(fc + i) ^ bj 
j =0 


k—2 

+ E 

3=0 


(h)s + f, g bj 


< b 2 VbVs 

~ |logfo| 2 | log b\' 


(3.38) 


Remark 3.9. The constants in Lemma 13.81 are independent of D,K , see also Re¬ 
mark 13.131 We need to keep track of the coupling between the modes in (I3.38p in 
order to study the linearized system close to bf and close the shooting argument, 
see (13.591) . 


Proof of Lemma This lemma follows from the orthogonality conditions (13.141) 
and the boundary conditions of lemma 13.51 


step 1 Computation of Mod. The Mod estimate follows from the sharp choice of 
orthogonality conditions (13.141) . Indeed, for any 0 < j < k, we take the scalar 
product of (I3.35|) with rj b j and use the orthogonality condition (I3.14|) to compute: 

-(e,d s ri bjj ) b + y(e, \y\ 2 Vb,j)b = (^^bj) - (a - b)(Ae, r] bJ ). 

We now intergrate by parts and use (13.141) again to compute: 


-(e,<9 sVb,j)b + \y\ 2 Vbj)b + (a - b)(Ae,ri b j) 


, * 

= 


= (e, —~ffbd b r) bj j + y \y\ 2 y b j + (a - b)[-Arj bJ + by 2 rj b>j ]) 

,2 


jo , by 
-b*%Vb,j + -2~Vb,j 


+ (a - b)[2bd b rj bJ - Arj bJ }) 


We evaluate all terms in the above expression. From (12.791) . (12.801) . (lA.IOp : 



to by 2 

-bo b rj b j + — 


<h| | log 6| ^ |<h|| log b\ V£ 

T l|e||b ^7T~ ^ 7T 
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Similarily from (12.82ft , (|2.83p . and (13.301) : 

|(e, (a - b)[2bd b ri bd - A rj btj ] 


< 


| b — a\ 


Vb 


Mlb ^ 


Vs 

Vb' 


We now estimate the T terms given by (13.361) . From (13.341) . (12.831) : 
\{^,Vb,j)b\ ^ ll^llblkbjllb ^ 


6§ 


$1 


|log6| V b 


|log&| <h |4>||log 


Vb 


The collection of above bounds together with (12.831) and (13.361) yields 


| (Mod, 77bj )b| 


< 


b , l$l|logfr| 

\ + 1^11 log 6 |‘ 

Vs 

6 

[ 6 2 J 

Vb_ 


(y b ,j,Vb,j)b ~ |iog6| 2 
We now argue differently depending on k. 

step 2 Case k = 0. In this case, b = 6q and thus from (13.26|) : 


4* — (Po)s + 26o(a — bo) — b s + O 


+ 


VbVE \ 


log b | 2 | log b | ) 


This together with (13.25}) implies the bound: 


|<F| < 


log b| 


+ \b 8 


and thus ()3.39|) . (I3.32p . and (I3.26P imply: 


I b s + b 2 X bt o\ 


< 


b 2 + I bs 


|log6| 2 |logb| 
Using the rough bound (|3.25p . this gives: 

2 b 2 


+ 


1 + 


\b s \\ log6| 

b 2 


l + O 


|log6| 


+ 


|logb| 


l + O 


|log6| 


VbVS 

|log6| 2 ’ 

~ |logi-l 2 


and (13.371) follows. 


step 3 Case k > 1. In this case, we have from (|3.9D : 

\b s \<b 2 

and we therefore need to estimate <h: 

k 

= b s + 2 b(a — b) = b s + 2b 


+ 2b'^2(bj - bj)( 1 + 


3 =0 
k 

2b ^2'bj + O 

3=0 


2a j 
|log6| 


2 a. 


j=o 


+ 26 


1 + 


|log6| 


-6 


b 2 


VWs\ 

|log6| 2 ' | log 6| J 


+ 


where we used (I3.26p . (13.121) and the bootstrap bounds (13.231) . (|4.13D in 
step. This implies in particular the rough bound 

6 2 


m < 


log 6| 


(3.39) 

(3.40) 

(3.41) 

(3.42) 


(3.43) 
the last 


(3.44) 
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From (13.391) it follows that: 


| (Mod, ?7b,j)bl < b 2 VbVS < ft 2 


(Vb,j,Vb,j)b ~ |log6| 2 |log6| 2 ~ | log b\ 2 ' 

We now recall (13.1 ID and compute from (13.261) . (13.121) : 


(3.45) 


a - y ~](b) + bj) 

1=0 


1 + 


— n e + ^ ] bj + O 

i=o 


2a l 

log b| _ 

6 


+ 0 


+ 


+ 




\J log 6| 2 | log 6| a/6 ) 

Vs \ 


I log £>| 2 \\ogb\Vb ) 
We now use (13.431) . (13.12ft . (13.461) to compute explicitly: 


(3.46) 


, , , , 2 (a — b)b k kb k ^ 

(b k ) s + bb k \ b)k + | bg6 | + —^ 


— {b% + b k ) s + b(b% + b k ) 
2(a - a e )b k 


2 k + 


log b\ 


+ o 


1 


+- 


log b\ 

1 


+ k(b1 + b k ) 


2 yz~ b j+° 

3=0 


log b\ 2 
b 


+ 


2 (a*-b){b% + b k ) 


log b\ 


+ 


Vs 


|log6| 2 Vb\\ogb\ 


— ( b k )s H— 

s 


b k + (k + 1) bj 

j=0 

Similarity for j = k — 1: 


+ o 


b 2 VbVS 


|log6| 2 | log 




|log6| 


— (b k —i)s 4~ bb k —i 


2(a - ci e )b k —i 


2(k — 1) + 


log 6| 


+ 0 


log 61" 


+ 


2(a e - b)b k - 1 


|logb| 


|log6| 


— (b k — i)s + 

s 


(k - l)b k -\ - k(b% + b k ) 


2 bj +o 

1=0 


b + VS 


|log5| 2 \/6| log 6| 


lfe-1 


5 fc _i - (fc + 1) 

l=o 


+ o 


b 2 VbVs 


\\ogb\ 2 |log6| 


Finally for 0 < j < k — 2: 


(6,)» + bbj\ b j + tytyty + A-0 + l)^i „ 


J 

(bj) s + + 0 ^ji 0 gk|2 1 |i 0 g5| J 


|log6| 

b 2 


+ 


VbVs \ 


Injecting the above bounds into (13.451) yields (13.381) . 


□ 


3.4. Energy bound. We now arrive at the second main feature of the analysis 
which is the derivation of suitable energy bounds for e. The key here is the dissi¬ 
pation embedded in the problem and its geometry which feeds back into the energy 
estimates through the boundary conditions. A careful analysis of this interaction 
will allow us to close the energy estimates. 
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Proposition 3.10 (Energy bound). There holds the pointwise control: 
1. for k = 0: 


1 d_ 

2 ds 


\s + o 



+ cbS< 


6 4 

|log6| 2 ’ 


2. for k > 1: 


1 d_ 

2 ds 


| S + O 


|log6| 5 / 4 ) } 


+ [3k + c\b£ < 


I<b A 

|log6| 


for some universal constant c > 0. 


(3.47) 


(3.48) 


Remark 3.11. The sharp coercivity constant 3k in (I3.48|) which follows from the 
sharp Poincare estimate (1A.3I) is essential to close the energy bound, see (13.5811 . 

Proof of Proposition \S.1CH We compute the energy identity for £ and estimate all 
terms. 


step 1 Algebraic energy identity. Recalling from (13.611 that £2 = %&£, it follows 
from (13.351) : 

d s S2 + TL a £ 2 = [d s , 7 ib}s + [H a , / Hb\£ + HbJ 7 - 
To compute the commutators [d s ,'Hb[, [Ha^b] we use 

[A, A] = 2A (3.49) 

which yields: 

[d s ,7i b \£ + [H a ,'Hb]£ = b s Ae + [Hb + (a - 6)A,'H b ]e = b s Ae + (a - b)[ A, - A] 
= b s Ae + 2(a — b) Ae = (b s + 2 6(a — b)) Ae — 2(a — b )[—Ae + 6Ae] 

= <f>Ae + 2(6 — a)e2- 

Hence the £2 equation: 

d s £2 + TLa £2 = 4>Ae + 2(6 — a)s 2 + TLb^~ ■ (3.50) 


We now compute the modified energy identity: 


1 d _ Id 

2 ds 2 ds J y >\ 2 


j. 


e\e 2 ydy = —- 


by 


b\v\ 


y |e 2 | e 2 ydy + (<9 s e 2 ,£2)& 


'2/>l 


— —^ || 2/^2 life + (4>Ae + 2(6 — a)e 2 + Tib^~ ~ kL a £ 2 , £2 )&• 

We carefully integrate by parts to compute: 

.,2 


'?/>! 


by^ C by ^ 

£ 2 pi a £ 2 e 2 -ydy = - £ 2 [H b £2 + («• - 6)Ae 2 ]e 2 “ ydy 

Jy> 1 


= [ dy(pbydy£ 2 )£ 2 dy + (6 - a) [ e 2 y9 y e 2 e 2 ydy 

Jy>l J 

f n by 2 

= -p b {l)£2{l)dy£ 2 {l) - / I^£21 e 2 ydy 

Jy>l 


b CL /-, \ 9 /i \ ^ ^ 

-p fe (l)ei(l) - —- 


[ 4 [ 2 - ^2/ 2 ] e ^ ydy 

Jy> 1 


12 1 / 7 \ 11 112 11 112 


= —IIV 2 IIS + (a- 6)||e 2 ||i - 


1/3! ||& - Pfe(l)e 2 (l) 


o b - a 

9y£2 + ^ £2 


(I)- 
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This yields the algebraic energy identity: 

= -\\d y £2\\ 2 b ~ (a - b)\\£ 2 \\l - ^\\y£2\\ 2 b ~ Pb£2 


\ b-a 
9 y £2 + ^—£2 


+ ^(Ae, £2)5 + {T~Lbd~ j £2 )b- 

We now estimate all terms in the right-hand side of (13.5111 . 
step 2 Nonlinear estimates. From (I A. 101) . (I3.44|) . (I3.25|) : 

\$\\(Ae,e 2 ) b \<J^^V£ < b 7 £ 


|logb| b ' ~ ~ |log 6 |' 
Moreover from (13.1411 (T^Mod, £ 2 )b = 0, and from (13.3411 . (13.4411 : 


I (n b *,e 2 ) b <bV£- 


bl 


|log 6 |' 


We now estimate from (|A.12(1 . (13.4411 . (13.27f) . (13.3011 : 

b 2 


\nye 2 Wl 


\\dyS2\\l 


+ b 6 


\\dye 2 Wl 6 4 


< w 


+ 


|log 6 | |log 6| 2 


(1) 

(3.51) 


|log 6 | [ b 2 

It now remains to treat the boundary term in (13.5111 and we argue differently de¬ 
pending on k. 


step 3 Conclusion for k = 0. We compute from (13.2811 . (13.371) : 


d y £ 2 ( 1 ) = —a s - 


2 b 2 
|logfo| 


+ 0 


= -(a-b) 3 + 0 
and hence using (13.271) : 
p 6 (l)£ 2 (l) 


o b-a 
dy £2 + ^—£2 


( b 2 
V Uog b\‘ 
f Vby/£ b 2 

\ |log 6 | + |log 6| 2 

(1) 


= -(a - b) s - (b s + 


2 b 2 
|log 6 | 


) + o 


b 2 


|log 6| 5 


= e 2 a (a — b) 


a(a-b) 2 

— (a — b) s -|-—-1- O 


( Vby/£ b 2 


y |log 6 | |log 6 | s 


= e 2 [~( a - b) 2 (a - b) s - b(a - b)(a - b) s ] +0 (b ^ ^ + 


|log 6 | |log 6| 2 


d f _b(a — b) 3 ^ _bb(a — b) 2 \ b s e 2 (a — 6) 3 


ds 


6 


b b(a — b) 2 
+e 2 ’ 


4 . — 

2 b 


+ 0 b 


2 6 
+ 


bz\f£ 


|log&| |log 6| 2 


d f _b(a — b ) 3 _bb{a — b ) 


ds 


e 2 


6 


+ e 2 . 


+ 0 b 


b^VS 6 4 b 2 (a — b) r - 
|log 6 | + |log 6| 2 + |log 6 | 


We now observe from (I3.26p . (13.251) that 

\a~ b\ < 


|log 6 | 


(3.52) 
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and hence the collection of above bounds yields the control: 


is 1 -'*" 




b^VS 6 4 


|log 6 | |log 6| 2 j ' 


We now inject (13.271) . (1A.11D with k = 0 and (|3.47|) follows, 
step 4 Conclusion for k > 1. Let 

k 

a = a e + ^2 bj 

j =o 

be the leading order part of a, see (I3.46|) . Then from (13.381) . (I4.13|) . (13.121) : 

k 

's 


1 + 


2 a 


■j 


(is + 'y ' ^b,jbb 3 i .. ,, 

k k 

= {a e ) s + ^2(bj)s + X/ ^ b ’jb(bj + bj ) 


= | a e - 6 £ 1 + 


3 =0 1=0 

“ 2 otk 


i 2a,- 
1+ J 


|logfo| 


+ O 


I log 6 | 

/ Kb 2 vWf\ 
y |log 6| 3 / 2 + |log 6 | J 


+ 1 + 


= -2 


_ 2 afc\ 

I log b\ J 

(a e - b)b k e 

I log b\ 

O ( Kb ' 2 ^bV£\ 

2 fcs 2 |logs| \Jlog &| 3 / 2 |log 6 | J 

We therefore estimate using Lemma [3.51 


[(bt)s + \jbbl\ + O 

( Kh 

e +o 


b 2 


I log b\ 2 

( Kb 2 VWe\ 

1 |log 6| 3 / 2 + |log 6 | J 


k + 1 


~Pb£ 2 


o , b-a 
d y £2 + ^ £2 


(1) 


= —e 2 a(a — b) 
= —e _ 2 a(a — b) 


a s T y ( \fyjbbj 

3=0 


i 2a,- 
1 +- - 


|log 6 | 


+ o 


(— 

v |iog&h 


f Kb 2 VWe\ 

2ks 2 \ogs ' \Jlog 6| 3 / 2 |log 6 | J 


( ~\ k + 1 ^ 

(a - ajs - 01 - 1 - O 


= —e 2 a(a — b)(a — a) s + e za(a-b)- 


k + 1 | c ( Kb 4 | bb?Vs 


+ 


2ks 2 \ogs l |log 6| 3 / 2 |log 6 | 
= —e _ 2 (a — a) s [(a — a ) 2 + (a — a)(2a — b) + d(a — 6 )] 


+o( 

b ) 
|log 6 |' 

d 

1 e- 

ds 


+°\ 

( b 4 
l log 6 


Kb 4 662 Vf 


y|log 6| 3 / 2 |log 6 | 
(a — a ) 3 (a — a ) 2 


+ 


( 2 d — 6 ) + (a — a)a(a — b) 


+ 


Kb 4 ,&&§>/£, 3 


+ 


|log 6 | |log 6| 3 / 2 |log 6 | 


+ Kla — a I 
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We have the bound from (13.461) : 

a — a = O 


+ 


Vs 


|log6| \logb\Vb I 


Injecting the above bounds into (13.511) and using the rough bound (|3.25|) we obtain 
the bound: 


1 d 

2 ds 


S + O 


b 3 


|log&| 


= -|M2" (a -&)N||I 

(\\dyS2\\l + |N|b , ,b^VS 


+ o 




|log6| 


^ + b-— + 


b A 


|logb| |log6| 


We now estimate from (13.461) : 


a — b = -b O 

2s 


|logfo| 


= kb + O 


|log6| 


and (13.481) now follows from (1 A. Ill) with (13.271) . 


(3.53) 

□ 


3.5. Proof of Proposition 13.21 We are now in position to give a sharp descrip¬ 
tion of the singularity formation for our set of initial data. The key is to close the 
bootstrap bounds of Proposition 13.21 We distinguish the cases k = 0 and k > 1. 

step 1 Closing the bootstrap bounds. Our goal is to show that the bounds (13.211) . 
(13.221) improve in the case k = 0 and similarly for the bounds (I3.2ip . (13.231) . (I4.13|) 
in the case k > 1. The improvement of the energy bound (13.211) will follow from 
proposition 13.101 while the bounds (I3.23p . ([4.131) will be improved for a suitable set 
of initial data constructed via a topological argument. 

k = 0 . First observe that (13.371) ensures 

b s < 0 and hence b(s ) < bo < b*. 

From (I3.52p . 


As 

A 


+ b = b — a = O 


and hence: 


logA(s) = - [ b 
Jo 


1 + 0 


|log6| 


|log6| 

da < +oo implies A(s) > 0. 


(3.54) 

(3.55) 


We now rewrite ()3.48l) (with k = 0) as 


ille + o 


ds I 2 


b 3 


|k>gfr|- 


+ be 


S + O 


b 3 


with c > 0. Using (|3.54p we obtain: 


d f 1 


ds I A c 


1 


-S + O 


b 3 


|log6| 2 

< ^ 

~ A c |log6| 2 ’ 


~ |log6| 2 


(3.56) 


(log^l 2 . 

We now integrate in time. To evaluate the right hand side, we integrate by parts 
using (13.371) : 


l 


b A 


o A c |log6| 
1 b 3 
c A c |log6| 2 


;da = 


/'[ 


A. s 


b 3 


A A c |log6| 2 


i rbs 

c Jo A c 


3 b 2 


+ 


+ o 


2 b 2 


b A 


llogfe^ 1 


da 


|log6| 2 |log6| 3 


da + O 

Jo 


A c |log6| 5 


da. 
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Using the smallness of 6 we get: 


f 


r da < 


b 3 


r(»)- 


/ 0 A c |log 6| 2 ~ A c |log 6| 2 
Hence (13.ID . (13.151) and the time integration of (I3.56|) ensure: 
b 3 


~ 11_LI2 ( S ) + 


|log 6 | 

We moreover estimate from (13.371) : 
d 


h 3 

£( 0 ) + 0 


|log 6 0 | 2 J ~ |log 6| 2 


<-^W + A'W b ° 


|log 6 0 | 2 ' 


(3.57) 


1 <- 3 1 

6 3 r 

i A c |log 6| 2 J 

A c |log 6| 2 


36, 26 

+ 


A, 


— c- 


6 6 |log 6 | A 


6 3 


A c |log 6 | 


c 6 + 0 


|log 6 |- 


> 0 


and hence using (13.11) again: 


A c (s)rri < A c (so)t 


|log 6 0 | 5 


|log 6 | s 


which together with (|3.55l) implies 6 (s) > 0 and closes the bound (13.221) . Injecting 
this into (13.571) improves the energy bound (13.211) for D universal large enough, 
which concludes the proof of Proposition 13.21 for k = 0. 

k > 1. This case requires a shooting argument to build the nonlinear manifold of 
perturbations (I / j)o<j<fc- 1 ■ We first rewrite (|3.48j) using (|3.9I) : 

6 3 


d_ 

ds 


£ + 0 


Using (13.151) . an integration-in-time yields: 
3+1 


£(s) < 


s 3 +t 


£o 3- 


|log 6 | 
e 

6n 


+ 3 + |j l -£<^— 
ki s s 4 1 logs| 

s Ko 3+ £ 


K 


< 

~ s 3 (logs) 


< K 


log 6 0 | 

6 3 


6 3 1 

+ 77-77 + 


log 6 | s 3 +f J s cr 4 |logcr| 


/’ 

-'SO 




|log 6 |' 


(3.58) 


This means that there exists aC>0 universal large enough such that if D = CK, 
the bound (13.211) gets improved, and we assume it now. We inject this relation into 
(13.381) and conclude: 


{bk)s 3 — 
s 


bk + (k + 1 ) bj 
3=0 




~ 1 


+ 

(Ofc-i)s 3— 
s 


k~ 1 ; 


J k -1 


{k + l)Y^ bj 

3=0 


k—2 

+ E 

3=0 


(bj)s + 


<vm^< 


Vk 


|log 6| 3 / 2 ~ s 2 (logs ) 3 / 2 
Equivalently using the change of variables (13.161) : 



+ 

(Vfc-i) s -|— 

k 





1=0 



k—2 


+ E 

j=0 






(3.59) 
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The bootstrap bound (14.131) implies that \Vj\ < 5K, j = 0,... k — 2. Therefore, 
from the first two bounds in ()3.59p we conclude that for a sufficiently large K the 
following bound holds 

SK 


(V k ) s + ^Y(V k + V k _ i) 


+ 


(V'fc-i)s- (V k + (1 + d k )V k _ i) 


< 

r^j 


where we remind the reader that d k = jTfrppj' Recalling the definition (13.171) of the 
matrix A k , the above inequalities can be succinctly rewritten in the form 

' 5K \ 


d. 


V k 

Vfc-i 


k+ 1 


Ai. 


V k 

V k -1 


+ o 


s J 


which in turn leads to 


(w^)* + {k + 1) ^ w k 


+ 


(Wfc_ i) s + (fc + 1) ^ lTfc-r 




(3.60) 


where W k _\ are defined in (13.181) and < 0 < are the eigenvalues of A k . 
This first yields the control of the stable direction W k , since after integrating-in-time 
the first bound in (I3.60p . we arrive at 


S k + i)mJ 


\W k (s)\ < \W k (0)\ 


+ 


s {k+l)fi% s (k+l)fi 


i r 

^Js 0 


s < 57TcrC fc + 1 )^i 


(7 


dcr < 1 + C^iL, 


where we used (I3.19P and the positivity of /jl\. This improves (13.231) for K suffi¬ 
ciently large and 6 < We now argue by contradiction and assume that for all 
(■y-,..., Vk y 2 , W k -\) € the boostrap time s* is finite, so that from (14.131) : 


k—2 


iWk^is*)] 2 + J2 

1=0 




= K 2 . 


(3.61) 


We claim that this contradicts the Brouwer fixed point theorem. Indeed, us¬ 
ing (I3.59p . (13.601) . and the strict negativity of /Jy. 



|W fc _ 1 (s )| 2 + ^ 


k -2 T 7 - / \ 2 

Vj{s) 


1=0 


<5 

k—2 


(O 


\^\(k + mtys*) + Y ^vf{s*) + o ( SK 2 + ir 3 / 2 ) 


l=o 


for some universal constants c,C > 0. Hence 

k—2 


Vj(s*) 


(s*) > 0 


* i^)i 2+ E 

( 1=0 

for 0 < 5 < 1 universal small enough in (|4.13l) . Let 

T> ,V0 Vk~2 jjr \ 

V = j-,W k _ i), 

then this implies from standard from standard argument that the map 

Bx(R d ~ 1 ) 9 V(0) ^ s* 


(3.62) 
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is continuous, and hence the map 


Bk^ 1 ) - 
V'(O) eA V [s*(F( 0 )) 


-I'l 


is continuous and the identity on the boundary sphere 1 (-^ 0 , a contradiction 
to Brouwer’s fixed point theorem. This concludes the proof of Proposition 13.21 for 
k > 1 . 


Remark 3.12. Note that (Vo(To), • • • -I4_ 2 (^ 0 ), Wfc_i(eo)) are by construction lying 
on a nonlinear codimension k manifold of initial data. The fact that the set of such 
initial data forms a Lipschitz manifold in the H 2 topology reduces to a local unique¬ 
ness problem in the class of solutions satisfying the a priori bounds of Proposition 
13.21 Such a uniqueness problem has been recently solved in a related framework in 
the more complicated case of the wave equation [5] and the KdV equation [29], see 
also [?], and a completely analogous approach can be applied here. We therefore 
omit the details. 


Remark 3.13. Note that the presence of \[K on the right-hand side of (I3.59P 
is essential to the closure of the estimates. It originates from the bound (13.3811 . 
where we carefully tracked the constants and proved that only \f£ appears on the 
right-hand side of (I3.38[) . 

step 2 Global H 2 control. From Proposition 13.21 the solution remains in the boot¬ 
strap regime of Proposition 13.21 as long as it exists in H 2 which requires: Vs > 0, 


ll M ( s )llL 2 (|a;|>A(s)) + II Vlt(s) ||l 2 (| x |>A(s)) + ll^ w ( s ) IIl 2 (|x|>A(s)) < +°° (3.63) 

and 

A(s) > 0. (3.64) 

I 11 the case k = 0, the positivity of A follows from the time integration of (13.3711 
which implies 

|logA(s)l ~L b { 1+0( iif 0 da K +o °’ 

while in the case k > 1 we use a = —X s /X and the estimate (13. 5311 . which implies 
the above bound again. The global L 2 -bound follows from the basic dissipation law 
satisfied by the solutions of (EU): 

2^/"IHI i 2 (n(t)) + ll Vu lli 2 (n(t)) = °> 

which immediately implies that ||i*(s)||z, 2 (|a:|>A(s)) < 00 ■ The global iG-bound fol¬ 
lows directly from the dissipation of the Dirichlet energy m- For the global 
i7 2 -bound, we take a cut off function y = 0 for r < 1 and r = 1 for r > 2, then the 
weighted control (I3.2ip ensures 


H A *4 s )llL 2 (A(s)<r<2) < C(s) < +OO for S > 0 

since the exponential weight is uniformly bounded from below and above in A(s) < 
r < 2. To obtain the bound in the region r > 2 we compute: 


1 d_ 

2 dt 


J X \Au\ 2 = J x^dtuAu = j xA 2 uAu 


f X |VA 


ur + 


Ay|Au| 


(3.65) 
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and hence 


Jx\Au\\s) < || X Au(0)||| 2 + ^ 2 ^y||Aii(s)||| 2(A((T) < r < 2) dcT < +oo. 

Hence s* = +oo which concludes the proof of Proposition 13.21 


3.6. Proof of Theorem 11.11 We are now in position to conclude the proof of 
Theorem 11.11 


step 1 Finite time melting. We claim that the solution melts in finite time with the 
law (11.41) . (11.51) as a consequence of the time integration of the modulation equa¬ 
tions. 

case k = 0: From (13.371) . (13.211) . we obtain the following pointwise differential in¬ 
equality for b : 

i '- + Sf = °(lkw)' (3 ' 66) 

We now follow HUES] to derive the melting speed of A and sketch the proof for the 
sake of clarity. Multiplying (13.661) by we obtain 


Mog6 


= 2 + 0 


1 


b 2 \ |log6| 

The primitive of is — l2Sk — I and therefore 

* 11“ II. II. 


¥ + i = - 2s 

b b 


+ o(f ’W~b\ dT) 

Jo |logb| 


which implies 


logs < logs 




Hence: 


b = — 


1 + log& 
2 s 


1 + 0 


1 


1 


|logb| 


dr 


1 + log b 
2 s 


Taking the log yields 

log6 = — log2 — logs + log(—log&) + O 


1 + 0 


|log6| 


|log6| 


which reinjected into (13.67|) ensures: 
logs 


b = 


2 s 


1 + 0 /l0gl0g£\ 

V lo s s / 


Injecting this into (13.671) again yields: 


b = — 


-logs + loglogs - log2 + 1 + 0 


logs loglogs 
2 s 2s 

Recalling from (13.261) that 


+ 0 


1 + 0 


logs 


b --^ + °{+EW 


(3.67) 


, log6 = loglogs — log2 — logs + O ([ 

V lo S s / 


(3.68) 


(3.69) 


( 3 . 70 ) 
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we conclude that 


which gives 


logs loglogs . 1 . 
-dogA), = —-^^+0(-), 


—logA = ^-(logs ) 2 - ^logsloglogs + O(logs), 


which in turn gives 

—21og(A 2 ) = (logs ) 2 

This leads to 


\/—21og(A 2 ) = logs 


1 _ 2 ^log » +0 

logs 

loglogs 


1 


1 - 


+ o 


logs 

= logs — loglogs + 0 ( 1 ) 


logs 

1 

logs 


from which 


V-2i°g(A 2 ) = _£_ e O(i) 
logs 


and hence from (13.7011 . (13.691) : 

-AA t = -^ = b + 0 

A 


1 


slogs ) 2 s 


= l£i£ 0(1) = -V2|logA 2 |+0(l) 


This yields the pointwise ode: 

_ e V2|logA 2 |+0(l) ( ' A 2^ = 1 

which integration in time yields: 

A 2 (f) = (T - t ) e -V^(T-t)\+o(i) 

and (El is proved. 

case k > 1: We estimate from (13.211) . (13.81) : 

A s k +1 k + 1 ( 1 

A 2 ks 2k 2 slogs \s(logs) 3 / 2 

and hence there exists c* = c*(u o) such that 

A: + 1 k + 1 

—logA(s) = —jj— logs - -^-loglogs + c* + o s _> +00 ( 1) 


or equivalently: 


We conclude that 


fc+i 

A(s) = c(u 0 )(l + o(l)) ^ 1 ° S ^ 1 —, c{u 0 ) > 0 . 


S 2fc 


T = 


r+oo 

/ A 2 (s)d, 

Jo 


and 


T-t = 


s < +oo 


fc +1 


(3.71) 


(3.72) 


(3.73) 


(3.74) 


fc + 1 

t= r”( C 2 + »(l))S25#lda = (^ + „(l)) (l0es) ^ 

ds J S (J k 


1 

Sk 


This implies 


1 (T — t) k . 

_ — i- - -~fc+r( c + °( 1 )) 

|log(T-t)| t 


s 


(3.75) 
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which together with (|3.74|> yields the melting law: 


A (t) = ( c*(u 0 ) + 1)) 


, . fc +1 

C T-t ) — 

|log(T-i)|^’ 


this is (| 1.5 1) . 


step 2 Non concentration of the energy. Pick R > 0 and a cut-off function 

0 for x < R 


xr(x) 


for x > 2 R. 

Then for t sufficiently close to the melting time T: 


1 _d 

2 dt 

s 


Xi?|Vtz| 2 dx = / x_rV«-V dtudx = /XrVu-V Audx 


= - I Au [xrAu + Vxr • Vit] dx 

1 


= - j Xr\^r\ 2 + \ J 


' V “l 2 4r (v 


dx 


(3.76) 


cr 


dx 


and hence the uniform bound on the Dirichet energy ensures: 

f T 

\/R >0, / x R |Au| 2 dx < +oo. 

do 

Hence for all 0 < t < T — t, 

/ r rt+r 

Xi?|Vii(t + r) - Vu(t)| 2 dx = J xr J (d t Vu)(a,x)d, 

< t f Xij|<9tVn| 2 dx < r / Xi?|VAu| 2 dx < C(.R)t, 
do do 

where the last estimate follows by integrating-in-f equation (13.651) with x = XR and 
using (|3.76D . Hence for all R > 0, Vu(t, x) is a Cauchy sequence in L 2 (|x| > 2R) as 
t —>• T. We conclude from a simple diagonal extraction argument that there exists 
u* € id 1 (M 2 ) such that 

Vi? > 0, Vii(t) —)• Vu* in L 2 (|x| > 2R) as t —>• T. 

Moreover, the uniform bound on the Dirichlet energy C3D ensures 
X7u* € L 2 , Vn(t) — 1 Vu* weakly in L 2 as t —> T. 

Pick now 

A (t)B(t), B 2 (t)b(t) = i|log 6 (t)| for k = 0 
A (t)B(t), B 2 (t)a(t) = A(f)|loga(f)| for k > 1. 

Note that in both cases we have B(t) 2 > 1. Then from (13.761) : 

J XR(t)\^ u (r)\ 2 dx < J | Vn(r)| 2 dx + J %R{t)\ Au(r)| 2 dx 

and hence integrating over [f, T) and using ()3.77|) . (13.791) . (11.31) : 

/x« W |Vu(r)| 2 dx- /x,( t )|V.1 2 dx ~ + J t j r>m \ Au ( T )\ 2dxdr - 

If k = 0, we use (11.41) to estimate 

T-t 2b(t)(T-t) 1 


m = 


(3.77) 

(3.78) 

(3.79) 


R 2 (t) 


A 2 (t) |log&(i)| 


0 as t —>• T. 
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The above limit holds since by (|3.72l) 

2 b(t)(T-t) 1 ^ 6(i)eV 2 l log ( T -*)l ^ b (t)e 2 V\W\ < i 


< 


< 

rs_/ 


\ 2 {t) |log 6 (t)| ~ |log 6 (t)| ~ |log 6 (t)| ~ |logfe(t)| 

where the last bound follows from (|3.68p and (13.7111 . If k > 1, then 

T — t a(t)(T — t) 1 
R 2 (t) A 2 (t) c|loga(t)| 

The above limit holds since by (13.721) 


0, as t —>• T, 


0 as t —>• T. 


a(t)(T — t) 1 
\ 2 {t) |loga(i)| 


(T -% (T-t) 

< |log(T-t)|-F- _ < 


1 


(T-t) fc +! 

fc +1 


|logo(i)| ~ l lo ga(*)l 


0, as t —>• T, 


|log(T-t)|' 

where we used (13.73|> and (13.7511 . Letting t —> T, we conclude: 


We now claim that 


J l v u*| 2 = Hm y Xi?(t)l Vn ( r )! 2 ^- 

(1 -Xi?.(t))|Vri(r)| 2 dx = 0 


lim / 

^ T 4{M>A(t)} 


(3.80) 


(3.81) 


which together with (13.8011 . (13.7811 concludes the proof of (11.61) . Indeed, in the case 
k = 0, from (12.7811 . (1A.10D . (13.211) . and (I3.79|) we obtain: 




(1 - XR{t))\Vu(t)\ 2 dx < [ \Vu(t)\ 2 dx 

J X(t)<\x\<2R(t) 


'i<\y\<2B(t) 


\\7v(t,y)\ 2 dy 


< P 2 bB 2 (i) 


= & 2bB 2 {t) 


/ \bd y rib 0 \ 2 PbVdy + \\dys\\ 2 b 

Jl<y<2B(t) 


6 2 |log 6| 2 + 


6 |log 6| 5 


0 as t —>• T, 


|logb| _ 

and (I3.8ip is proved. A similar algebra holds for k > 1. This concludes the proof of 
Theorem 11.11 


4. Infinite time freezing regimes 

This section is devoted to the existence and stability of the freezing process 
emerging from strongly localized initial data. Throughout the section, we let 

± = +, P = P+, B> 0. 

4.1. Renormalized equations and initialization. We let 


u(t,x) =v(s,y), y=-—, A(0) = 1, 


(4.1) 


with the renormalized time 


s(t) = s 0 + f 
■Jo 


dr 

W) 


so > 1, 


(4.2) 
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and obtain the renormalized equation: 

d s v + 7iAV = 0, A = ^-, 
v(s, 1) = 0, d y v(s, 1) = —A. 

We now prepare our initial data in the following way. We let 

B ^ = 2? B&k = s fc+ 1 (logs)2 


so that with A e = B e : 


{Bl) s + BlB[2k + 2 + 


log s 


+ 


2(B-A°)B° k 
log s 


= O 


1 


s fc+ 2 (l 0 g s ) 3 


We define 

k 

Qp(y) ■= ~^2 B jVB,j(y) 

3=0 

and introduce the dynamical decomposition 

v(s,y) = Q/ 3 (s) +e(s,y), {f) B , j (s),e) B = 0 , 0 < j < k. 

We let again 

£2 = kL B £, £ := \\H B £\\ 2 b , 

which due to the orthogonality conditions (13.1411 is a coercive norm, see Appendix 
El We assume the initial smallness 

B{B%) 2 { 0) 


(4.3) 

(4.4) 

.(4.5) 

(4.6) 

(4.7) 


m< 


|log£(0)| 


and consider the bootstrap bound 


£ < 


DB(B e k ) 2 

| log B | 


(4.8) 


(4.9) 


for some large enough universal D = D(k ) to be chosen later. Moreover, we assume 
initially 

Bk(s o) = Bfao), s 0 > 1 (4.10) 

and bootstrap the bound 

|5 fc (s)| < 10 B%{s). (4.11) 

For k > 1, we also let 

Vj(s) 


Bj (s) = 


g/c+i (log s) 2 


j = 0 ,..., k - 1. 


and assume 


We define 


s = 


]T|^(s)| 2 <A' 2 . 
j=0 

_/ sup^.. 0 ( (l4.9p . (14. ll[) hold on [s 0 ,s]} for k = 0, 


(4.12) 


(4.13) 


sup„^„ 0 { (l4.9|) . (14.111) . (I4.13P hold on [sq,s]} for k> 1. 


The main ingredient of the proof of Theorem 11.21 is the following: 


Proposition 4.1 (Bootstrap estimates on B and e). The following statements hold: 

1 . Stable regime; for k = 0, s* = +oo. 

2. Unstable regime; for k > 1, there exist constants K, 5 = 5(K) <S 1 and 
(Vq(0 ), • ■ ■, 14 _ 2 ( 0 ), 14 _i( 0 )) depending on e(0) such that s* = +oo. 
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From now on and for the rest of this section, we study the flow in the boostrap 
regime s € [so,s*)- Note in particular the rough bounds 


\Bj\<\B e k \, 0 <j<k and £< 


V\ l °gB\ 


(4.14) 


for sq > s 0 (K) large enough. 


4.2. Extraction of the leading order ODE’s driving the freezing. We start 
with the constraint induced by the boundary conditions: 


Lemma 4.2 (Boundary conditions). There holds: 
k 


A ' liogfll 

e 2 (l ) = A(B-A), 


l + J%: + 0 


1 


I logS | 2 


+ o 


Vs \ 

\\ogB\VB) 


1 + 


<B e k . 


d y e 2 (l) = A S ~ A\B ~A)+Y, A B jBB 3 

3=0 

Remark 4.3. Note that (I4.15p . (14.141) imply 

^ ~ Bk + \\ogB\V~B 
Proof of Lemma 1,9.51 We compute from (12.891) (13.291) : 

d yVB,j( 1 ) = d y r]B,j(i)e~% = 1 + + o 

This implies that 

k k 

dyQfi{ 1 ) = — J2 B jdyVB,jV) = ~ Bj 

3=0 j =0 

Since v = Qr + e, it follows that 


2a i 


| logS | 


+ o 


| logS | 2 


| logS | 2 


(4.15) 

(4.16) 
<:.17) 

(4.18) 

(4.19) 


1 + |logS| + ° 


A. s 


| logS | 2 


-- y \y=l 


= V„ 


i d y Qp\ =1 — d y Qp{ 1) 


+ X! B j 

3=0 


'y =1 A 

2 a 

1 + 


| logS | 


+ o 


| logS | 2 


which together with (IB.21) yields (14.151) . From (14.31) . v(s, 1) = 0 and d y v(s, 1) = 
“X = ~ A: 

0 = H a v{ 1) = {U B v + (S - A)Au)(l) = e 2 (l) + A(A - B), 
this is (I4.16|) . Now from d y v(s, 1) = — A, we have 

d s d y v(s, 1 ) = -A s . 

On the other hand, taking d y of (14.31) . we have: 

0 = d s d y v + dy(H B v + (S - A)Au) = d s d y v + d y e 2 + dyTLsQp + (S - A)yAv. 

We evaluate the above identity at y = 1. From (I4.3[) and d s v(s, 1) = 0, d y v(s. 1) = 
-A: 


Au(l) = — AAv{\) = A 2 . 
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By construction, 


dyV-bQp — ^ ' ^B,j B Bjdyf/Bj , 

i =o 


and hence: 


- A S+ d y£ 2 (l)-Y^ XB,jBBj 

3=0 


1 + TT%-\+° 

| logS | 


logS| 


+ A 2 (B- A) = 0. 


□ 


We now compute the leading order modulation equations. 

Proposition 4.4 (Leading order modulation equations). Under the a priori bounds 
of Proposition \f.l\ there holds 

dsQp + T-LaQ/3 = Mod + 'L (4.20) 

where we defined the modulation vector 

k 


( Bj) s + BBjXbj + 


Mod := — 

3=0 

and the remaining error satisfies the bound: 

1 _ 1 

S 




^||s + —W'Hb'&Wb 


2 (S - A)Bj 
| logS | 

Vbbi 


VB,j, 


< 

~ | logS 


(4.21) 


(4.22) 


Proof of Proposition \3.(A Let the deviation 

4- := B s + 2S(S - A). 

By definition 

n A = n B + {B - a)a 

and we therefore compute from (14.6p : 

~dsQp( s ){y) ~ UaQp 


k r 

E 

3=0 

k 


B • 

(■ Bj) s fj B ,j + B s -ABd B fiB,j + BBjX B ,jfjB,j + (S - A)B j Afj B ,j 

JD 


y l [( Bj) s + BBjXbfifjBj + (S - A)Bj[Af)B,j - 2 Bd B f} B ,j] + ^-Bd B riB,j$ 


3=0 

k 

E 

3=0 


(Bj) s + BBjXb.j + 


2 (S - A)Bj 
I log S | 


VB,j 


+ (S - A)Bj ( A f) B j - 2Bd B fiB,j - BdsfjBj } ■ 


We now estimate from (14.1811 : 


and hence using (|4.14[l . ()2.93p : 


1*1 < \AB\ < BB\ 




* || < I lo gS 

\ T $Bd B r, B J B < -j=- 


logS| 


(4.23) 
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and similarily for higher derivatives. Moreover from (I2.95p . (|4.14|i : 


(.B - A)Bj 


- 2BdBVB,j — 


]VB,j 


< 


1 


B 


Vb I logs I 


BBl < 


y[BB i 


3 |logS| 

and similarily for higher derivatives, and (14.2211 is proved. 

4.3. Modulation equations. The relations (14.31) . (14.201) yield 

d s e + PLa £ = B. T = —Mod — 4/ 
and we now compute the exact modulation equations. 

Lemma 4.5 (Modulation equations for Bj). There holds for 0 < j < k: 

< bb% 


logS | 


□ 


(4.24) 


(Bj) s + BjB I 2j + 2 + 


I logS| 


logS| 


(4.25) 


Proof of Lemma M Let 0 < j < k and take the scalar product of (14.241) with f/Bj 
and use the orthogonality condition (14.7|) to compute: 

\y\ 2 

— ( £,d s fjB,j)B - B s (e, -^-fj B ,j)B = + (A- B)(Ae,f/B,j)- 

We integrate by parts unsing (HZJ): 

\y\ 2 

-(£,d s f)B,j)B - B s (e, -y- VB,j)B + (B — A)(A£,fj B ,j)B 

\v \ 2 

= — (e, d s fj B j )b ~ B s (e, - (B — A)(e, B\y\ 2 i] B ,j + Af) B ,j)B 

= -(£, l -(B s + 2 S(S - A))\ y \ 2 fj Bd ) B - (£, Bs + 2B £ B ~ A) Bd B fj B j)B 
+{B — A)(s, 2Bd B fiB,j — Af}B,j) B 

<J> 1 o 

= -(e, — \-B\y\ 2 fiB,j + Bd B f}B,j]) B + (S - A){e, 2 Bd B r)B,j ~ Ai] B j)b- 

JD Z 

We now estimate from (IB. 2D . (|2.92[) . (12.931) . (14.231) : 


-^,^[\ B \y\ 2 VB,j + Bd B f]B,j})B 


and using (12.951) and (14.7f) : 


,, I, |$| | logS, 

^ Mb -^ 


Vs 


I (s - A)(e,2Bd B f)B,j - A VB,j] 


b | < 


S y/B ~ | log B\VB 

Vs 


B 


\ £ \\b < 


Vb ' ~ | log s | Vb ’ 

We now estimate the T terms given by (|4.24p . From (|4.22D : 


\('&,VB,j)B\ < W^WbWvbjWb 


< 


VBBl\}ogB\ < 

| log B\ VS ~ k ' 


The collection of above bounds together with (12.961) and (I4.24|) yields 


| (Mod, 77bj)b| <• S 


(VB,j,VB,j)B 
or equivalently: 


Vs 


E 

3=0 


|io g s | 2 | log s| Vb 
2 (S - A)B 


+ Bi 


BBl 


logS| 


( Bj) s + BBj \ B) j + 


logS| 


< bb% 


log S | 2 


(4.26) 
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We conclude from f|4.26|> . (|4.18p . (I2.89f) : 

4 


( Bj) s + BjB 

this is (|4.25p . 


2j + 2 + 


log B 


+ o 


log B | 2 


. BB% 

~ | log ,B | 2 ' 


□ 


4.4. Energy bound. We now derive the energy estimate in the freezing regime 
Proposition 4.6 (Energy bound for freezing). There holds the pointwise control 
1 d 


2 ds 


£ + 0 


B I 


V I lo S 


e ^2 


+ B(2k + 4 + c)\\s 2 \\ 2 B <B^^ 


(4.27) 


for some universal constant c > 0. 


Proof of Proposition \4-0\ We compute the energy identity for £ and estimate all 
terms. 


step 1 Algebraic energy identity. Recall e 2 = BbE- We compute from (14.241) : 

d s £2 + TLa £2 = [d s ,H B \£ 4 - [Ha, HIb\£ + H. B 4F. 
and hence using (|3.49p : 

[d a , B b }e + [H a , U b \e = -B s Ae + \U B + (B - A) A, U B }e 
= -B s Ae + (B - A) [A, - A] = -B s Ae + 2 {B - A) Ae 
= -{Bs + 2 B(B - A))As - 2 {B - A)[-Ae - BAe] = -<4>Ae - 2{B - A)e 2 . 
Hence the £2 equation: 


d s £2 + HLae 2 — —4>Ae — 2 (B — A4)e 2 + PL b B ■ 
We now compute the modified energy identity: 


(4.28) 


-£ = 

2 ds 2 ds 

B K 


'y>i 


2 y 1 
E 2 e 2 y dy = 


ly>l 


y 2 M 2 e 2 ydy + {d s £ 2 ,£2)B 


\y£2 1| b + (~ 4*A£ — 2 (B — A)e 2 + HbB — PLa£2, £2 )b- 


We carefully integrate by parts to compute: 


/ • 
Jy> 1 


Byj 


j • 

Jy> 1 


By 


E 2 B.AE 2 E 2 ydy = - / s 2 [Hb £2 + (B - A)As 2 ]e 2 ydy 


By ; 


' 2 />l 


dy{p B yd y £ 2 )E 2 dy - (B - A) / £ 2 yd y £ 2 e 2 ydy 


-pB(X)£ 2 {l)dy£ 2 (l) ~ [ \dy£ 2 \ 

Jv> 1 


' J/>1 

O Byf , 
“e 2 yay 


' 2 />l 

+^-^Ps(l)e|(l) + B A I e\ [2 + 5y 2 ] e^ 2 “ ydy 
1 A Jy >1 

= — II^ 2 / £ 2 ||b + (H — A)||e 2 ||| 4 - 1 | 2 /^ 21 |s ~ / 5 b( 1 ) £ 2 ( 1 ) 


, B-A 

dy£ 2 -7i- £ 2 


(I)- 


This yields the algebraic energy identity: 

2~ds^ = -W d y £ i\\ 2 B ~{B~ A)\\e 2 \\l + j||ye 2 ||| - Pb£2 
— 4>(Ae, £ 2 )b + (' H b J r , e 2 ) B - 



B-A 


-£2 


( 1 ) 


(4.29) 


2 
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We now estimate all terms in the right-hand side of (|4.29l) . 
step 2 Nonlinear estimates. From (1B.12() . (I4.23J) : 

\<S>\\(A£,e 2 ) B \ < BB e k ^- < B- £ 


|lo g 5|’ 

Moreover from (|4.7[> . (14.211) . (T^Mod, e 2 )b = 0, and from (14.221) . (14. 14}) : 




bstA < 


I ) 2 


log B\ ~ | log 51 

We now estimate from (IB.141) . (I4.23[) . (14.16(1 . (I4.18|) : 


|d>||ye 2 ||| < BB% 


ll^ 2 ||| + £ |(i) 


B 2 


B 


\\dyS2\l 2 B + &A 2 


< 

rv .J 


II Va lip B 3 


+ 


( Bt ) 2 + 


| log B | | log B 

\\dy£2\\ 2 B+b B(B^ 


logB| |logB| 

£ 


B|logB| 2 


~ | log B 


log B 


2 • 


step 3 Boundary term and conclusion. It now remains to treat the boundary term 
in (14.291) . First from (14.161) . (14.171) . (14.181) : 


\B-A\\e 2 {l)\ 2 <\B\{AB) 2 <B\Btf 


Let 


A ~ Bj ( ! + 

3=0 


2a i 


logB|y 

and observe the bounds from (|4.15|) . (14.141) . (14.251) : 


\A-A\ 


< 


Bl 


log B | 


\As\ < BB% 


(4.30) 


We now rewrite (|4.171) using (|4.25p . (14.181) . (14.261) : 

k 


dy £ 2( 1 ) — (A —A) + £ 1 + 


j=o 


2a ,■ 


I log 51 


( Bj) s + A B jBBj 


+ o 


BB% 

log B\ 2 


= (M-A) s + 0 

from which: 


BB e k 

log 51 


p B (l)e 2 (l) d y £ 2 - 
= Pb (1)A(B-A) 

= o(b 


B-A 


-£2 


( 1 ) 


(. A-A) s + 0 


BB e k d2d6 
k + B 2 Bi 


(B(Bl) 2 


\ | log B\ 


log B 

+ PB (l )A{B-A){A-A) 
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We now compute using (14.301) : 

- Pb (1)A(B -A)(A- A)s = ps(l)(A - A) a [(A - A) 2 + (A — A)(A — B) — A(B - A) 




(A-A) 3 ( A-A) 2 (A-B) ~ ~ 

v ; + -- J -+ -- - {A - A)redA(B - A) 


~t" b( d 


3 2 

(A-A) 3 (A-A) 2 (A-B) 


+Pb(1) 


3 2 

(A - A) 2 


-(A- A)A(B - A) 


d_ 

ds 


°(S? 


fc\2 
e / 


SI 2 


(A s - B.) - (A - A)(A a (B - A) + (B s - A S )A) 

B 2 (B\ 


+ o 


n? 


log B | 


Injecting the collection of above bounds into (|4.29D yields: 

(B(B 2 ) 2 \ 


ki< e+ ° 


VI iog- 8 ! 7 


— —\\dy £ 2\\B — (B — A)\\£2\\ B + O 


(\\dyeg^ + B B(Bl) 2 


log B\ 


log B | 


and hence using the coercivity (IB. 13|) , (14.161) , (14.41) and (|4.18|> : 


1 d 

2 ds 


8 + 0 


B(B 2 ) 2 
I logB| 


+ B 


2k+ 5 + 0 


logB| 


l £ 2|| b 


2 < B(B t 


e\2 


I logB| ’ 


and (14.271) is proved. 


□ 


4.5. Proof of Proposition 14. ll and Theorem 11.21 We may now close the boot¬ 
strap estimates of Proposition 14.11 

step 1 Closing the energy bound. First observe = |A| < B and thus | logA(s) | < 
C(s) implies A (s) > 0 on [0, s*]. The control of the 8 norm easily implies the H 2 
control ||tt(s )||#2 < C(s ) and hence the solution is well defined from the point of 
view of the H 2 Cauchy theory on [0, s*]. We now integrate in time the bound (I4.27|) 
and obtain for some c > 0 

£ (.) < (?p i|0) + C «(f + C ifM»», 


log s\ 


g2k-\-3+c 


so 


logB| 


* '?) 


2fc+3+c 


8(0) + C 


B(s)(B*(s)) 2 . B(s)(B* k (s)) 


| logs) ~ I log sI 

Hereby we used the explicit formulas (|4.4|> to infer that s2 fe+ 3 + e Jq B \\ogB\ a 2k+3+c da < 
B(s)(B e (s)) 2 

- 1 log fc s | - and in the last inequality the initial data assumption (14.81) . This closes 

the energy bound (14.91) . 


step 2 Control of Bj. We estimate from (|4.25|> : 


from which 


{Bk)s + BB}. I 2k + 2 + 


d_ 

ds 


log B | 


< m 


log B\ 


(s k+1 (log sj z B k ) 


< - 1 -. 

~ s(logs) 2 
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An integration-in-time yields: 

Bk{s) = 


Bk{so ) 


+ o 


(4.31) 


s fc+1 (logs) 2 [B-(so) V lo g s o. 
and the initial data assumption (I4.1UI) now improves (14.111) . 

For k > 1, we now argue by contradiction and assume that for all ( Vj)o<j<k-i with 
Y^lj =o I^-(0)1 2 < K 2 , there holds s* < +oo i.e. 

k -1 

^|^( S *)| 2 = iF 2 . 
l=o 

We estimate using the variables (|4.12p and (14.251) : 


(ft), ~ 


Hence at the exit time: 

1 d 

2 ds 


k-j 


< 






k -1 

Era 2 

i=o 


ra)>^EiL(»*)i 2 >o 


l=o 


and a contradiction follows as in the melting case using Brouwer fixed point theorem. 
This concludes the proof of Proposition 14.11 

4.6. Proof of Theorem I1.2L The proof of Theorem 11.21 now follows from a simple 
time integration of the modulation equations. 

We estimate from (|4.25|> 


( Bk) s 4-- I k + 1 + 


log s 


B e 

< 


s(logs) 


and hence 


d_ 

ds 


{s k+1 {log sfB k ) 


< -I- 

~ s(logs) 2 


which implies for s large enough: 


_ c(u 0 )(l + o(l)) 
k s fc+1 (logs) 2 

for some universal constant c = c(uq). We conclude from (14.151) : 

A,, 4 c(u 0 )(l + o(l)) 


f"- 4 " 


q/c+I 


(logs) 5 


from which there exists Aqo > Aoo(rto) > 0 with 

( c(uo)(l+O a _> -|-oo M 

AqO A(s) — 1 


if k = 0, 


logs 

c(hq)(1+Qs-> + oo (D) if 

k > 1. 

s k ( log s) 2 — 


Since for large s > 1 we have 

~c[~ — 5 m and (11.91) follow. Finally, recalling 


.By. 


that f]Bj = e 2 r]B,j we estimate: 


l v ( 7 ?s,i )| 2 = [ \d y r]B,j ~ Byrisjl 2 e By2 ydy= I \d z i/j B ,j ~ z^B,j \ 2 e z 

r|>l Jy> 1 Jz>VB 

= |logH| 2 [j + o(l)] / \Pj - zPj\ 2 e~ z2 zdz = c k [1 + o(l)] |logH| 2 


I z> 0 


z dz 
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for some universal constant c k > 0. Note that we used (|2.27p . Hence using (IB.21) : 


fc-i 


/ Mv ~ B k r , B , k )\ 2 < |log£| 2 ^T£ 2 + ||Ve||| 2 < (B e k ) 2 . 

J j =o 

Therefore, since the self-similar rescaling preserves the Dirichlet energy, it follows 
that 

[ |Vtt| 2 = [ |Vu| 2 = B 2 [ |Vr) s , fc | 2 + Omf) 

Jn{t) J\y\>X J\y\>l 

which yields (11.101) . To prove (11.71) . note that by integrating (11.11) and using the 
Stokes theorem, we arrive at the following conservation law: 


d_ 

dt 


f / u(t,x) dx - n\ 2 (t) ) = 0, 

\Jn(t) J 


(4.32) 


which holds as long as u € L l (Q(t)). To see this and evaluate lim^oo ||tt||z / 1 (n(t)) 
first observe that v satisfies 

1/2 




/ 

dy>l 


|u| ydy < ||u| 


B 


- B v 


e 2 ydy 


'y> i 


< ——I 

~ Vb 1 


V\\b- 


On the other hand by (14.71) it follows that 

k 


11 ^ v - ' n in ii nil B I log B I 1 || .. \/~B 

v\\b < ^ \Bj\\\r]B,j\\B + ||e||s < , p , 2 —7^—I- -^\\ 7 ~Lb£\\b, < 


3 =0 


log B\ 2 y/B B 


log 51 


where we used ()4.4D , (12.961) , (14.121) , (IB.121) , and (|4.9|) . The two previous inequalities 
lead to 

. 1 

PUHn) < | logg | ass ^ 00 . 

Therefore 

Il' w llz, 1 (f2(t)) = ~^ 0 as t —> 00 

since |A(f)| remains bounded. It follows in particular that the conservation law (|4.32p 
holds and formula (11.71) follows. 

This concludes the proof of Theorem 11.21 


Appendix A. Coercivity estimates in the melting case 

This appendix is devoted to the derivation of various coercivity estimates in the 
melting regime 

± = -, p = p-, b> 0 

which are used along the proof. We start with the standard compactness of the 
harmonic oscillator. 

Lemma A.l (Weighted L 2 estimate). Let u,d z u € L 2 (M 2 ). Then \/k > 0, 

jz 2k u 2 zpdz< k J ( d z u) 2 pzdz + J u 2 zpdz. (A.l) 

Proof of Lemma \A.l[ Indeed, we use d z p = —zp and integrate by parts to compute: 

J ( d z u — 5z k u) 2 pzdz = J ( d z u) 2 pzdz + 5 2 J z 2k u 2 pzdz — 25 J z k+1 ud z upzdz 

= j(d z u) 2 pzdz + 5 2 J z 2k u 2 pzdz — 5[z k+1 pu 2 ]^ 00 + 5 J u 2 [{k + l) k ~ 1 zp — z 2k p\dz 
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and hence for 0<<5 = <5(fc)<Cl small enough: 

J z 2k u 2 pzdz <k J {d z u) 2 pzdz + J u 2 pzdz + J u 2 z 2k ~ 2 pzdz 

and (lA.ip follows by induction on k > 1. □ 

We now claim the main coercivity property at the heart of the energy estimate. 

Lemma A.2 (Coercivity of H b ). Let k E N and 0 < b < b*(k) small enough. Let 
u S Hp(r > Vb) satisfy 


u(Vb) = 0, (u, ifb,j)b = 0, 0 < j < k. 

Then the following inequality holds: 

\\H b u \\ 2 L 2 > ||Au ||^2 + ||(1 + z)d z u \\ 2 L 2 (A.2) 

p,b p,b p,b 

+ \\{l + z)u\\ 2 l2 + b\\ogb\ 2 (d z u) 2 {Vb). 

p,b 

Moreover, there exists a constant > 0 such that 


and 


\\d z H b u \\ 2 L 2 > 

p,b 


2k+ 2 + 0 



\\H b u\\ 2 L 2 

p,b 


c k b 2 \H h u(Vb )\ 2 


\\zH b u \\ 2 L 2 < \\d z H b u \\ 2 L 2 + b\H b u(Vb)\ 2 . 

p,b p,b 


(A.3) 

(A-4) 


Proof of Lemma I A. SI This lemma is a simple perturbative consequence of the har¬ 
monic oscillator estimate 112.811 . (12.3211 . and a careful integration by parts to track 
the boundary term in (lA.2j) . 


step 1 Proof of (IA.2I1 . Pick a small constant 6 > 0, then from u(Vb) = 0, we may 
integrate by parts and compute: 

\\{H b - 6)u\\ 2 L 2 = {(H b - S)u, (H b - S)u) b = \\H b u \\ 2 L 2 - 25(H b u,u) b + 5 2 \\u\\ 2 l2 . 

p,b p,b p,b 

We may now use the spectral gap bound (I2.32p with (IA.1H and conclude that for 5 
small enough: 


\\H b u \\ 2 L 2 > 

p,b 


2 U IIl 2 p,6 + IK 1 + z ) u \\l 2 — 5)u\\ L 2 ■ 

r p,o p,b 


We integrate by parts using the general formula 

(d z u, v) b = - Vbu(Vb)v(Vb)p(Vb) - (u , d z v) b - (u , l -) b + (u, zv) b 


(A.5) 

(A.6) 


to compute: 

\\H b u \\ 2 L 2 = ||Au||^ 2 -2(u zz + —,zd z u) b + \\Au \\ 2 L 2 

p,b p,b Z P,b 

= \\ a u \\ 2 L 2 - ( d z (d z u ) 2 , z) b - 2 \\d z uf L 2 + \\Ku \\ 2 L 2 

p,b p,b p,b 

= || Au ||^ 2 + b(d z u) 2 (Vb)p(Vb) + (( d z u ) 2 ,1) 6 + ((d z u) 2 , l) b 

p,b 

- ||Au||^ 2 - 2\\d z u\\ 2 L 2 + \\Au\\ 2 L 2 

p,b p,b p,b 

= ||Au ||| 2 +b(d z u) 2 (Vb)p{Vb). (A.7) 

p,b 

On the other hand, integrating by parts: 

(H b u,\ogz) b = (u, 1 ) b + Vbd z u{Vb)\og{Vb) 
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and thus: 


b\d z u(Vb)\ 2 < 


|log6| 1 


\\H b u\\l 2 +\\u\\l 2 


P,b J 


<_I_ 

~ |log6| 2 


ll^ll 


2 

L 2 ’ 

p,b 


where we used (IA.5I) . Together with (IA.7D . f)A.5lh claim (1A.2D follows. 

step 2 Proof of (IA.3D . (IA.4p . Dehne the radially symmetric function 

, f H b u(Vb) for 0 < ^ < y/b 
” (2) = \ H„u(z) for 2 > Vb 


(A.8) 


and note that v € 0 ■ Consider 


(v,Pj) o 


W:=v ^ 0 ww° p ' 


Then from 


/ \d z w\ 2 e zzdz>(2k + 2) / |u>| 2 e 2 zdz > (2k + 2) 

J z> 0 J z> 0 

On the other hand, from (12.27|) : 


z>y/b 


w\e 2 zcP. 

(A.9) 


P- - 


|log6| 


Vv 


< 


H 2 , 

p,b 


|log6| 


from which for 0 < j < k: 

(v,Pj)o = H b u(Vb) 
2 


f z 2 2 

/ Pje 2 " zdz + (id 6 u, ^b,j)b + O 

0 <z<Vb |log6| 


/||Pyu|| L 2 6 ' 
\ |logb| 


|log6| 


(u, X b jiJj b j)b + O ( b\H h u{\fb)\ 




+ 


\\H h u\\ Lli 

|logb| 


= O 


lb\H b u(Vb)\ 


+ 


\\H h u\\ L *. 

|log6| 


where we used the orthogonality ( u,i/j b ,j)b = 0, 0 < j < k. Therefore 

\\v-w\\ H i b <b\H b u{Vb)\ + l^-. 

Injecting this into (1A.9[) yields (jA.3l) . We now apply flA.ll) to v and conclude from 

(l£ll), dX8j1): 

\\zH b u \\ 2 L2 < \\zv \\ 2 L 2 +b 2 \H b u{Vb )\ 2 <\\d z v\\ 2 L 2 + |M | 2 2 + b 2 \H b u(Vb )\ 2 

p,b p,0 p, 0 p, 0 

< \\d z H b u\\ 2 L 2 +\\H b u\\ 2 L 2 +b\H h u(Vb )\ 2 <\\d z Hu\\ 2 L 2 +b\H b u(Vb )\ 2 

p,b p,b p,b 

and (IA.4P is proved. □ 

We now renormalize Lemma I A. 2 1 by letting e{y) = u{y/by) which yields exactly: 

Lemma A.3 (Coercivity of T-L b ). Let k G N and 0 < b < b*(k) small enough. If 
£ € Hp b (y > 1) satisfies 

e(l) = 0, (e, rj b ,j)b = 0 for 0 < j < k, 

then 

\\U b £\\ 2 b > ||Ae|| 2 + 6||9 y e|| 2 + 6 2 ||Ae|| 2 (A.10) 

+ 6 2 ||(1 + Vby)e\\l + b\\ogb\ 2 {d y e) 2 (l). 
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Moreover, 


%n b ef L 2 > 

p,b 


2k+ 2 + 0 


1 


|log&| /_ 


b\\H b e\\l - c k b 2 \n b £{ l)| s 


and 


b\\yH h e\\l<-\\d y H b e\\ 2 b + \H b £(l)\ 2 . 


(A.ll) 

(A-12) 


Appendix B. Coercivity estimates in the freezing case 

We now consider the freezing regime 

± = +, p = p+, B> 0 

and the operator 

H b = —A — A, u(Vb) = 0 . 

We claim the analogue of Lemma IA.2I 

Lemma B.l (Coercivity of H B )• Let k € N and 0 < B < B*{k) small enough. Let 


B\y\ z 

•0 bj = *pB,je 2 

and u € H^(r > y/S) satisfy 

u(y/B) = 0, (u, = 0, 0 < j < k, 

then the following inequality holds: 


(B.l) 


\\H b u\ 


L 2 


> 


\\Au\\ 2 L2 + \\(1 + z)d z u\\ 2 L2 


(B.2) 


+ ||(1 + z)u \\ 2 L 2 + B\\ogB\ 2 (d z u) 2 (yfB). 


J P,B 


Moreover, there exists a constant c k > 0 such that 

1 


\\d z H b u\\U > 


T z — 
L p,B 


2k+ 4 + 0 


\\H b u\\ 2 l2 - c k B 2 \H b u(VB)\ 2 (B.3) 


J P,B 


and 


\\zH B u \\ L 2 < \\u z 

p,B 


\logb\J 

H b u\\ 2 L 2 +b\h bU {Vb)\ 2 


P,B 


(B.4) 


Proof. We follow the proof of Lemma (1A.3[) . 


step 1 Proof of (IA.2I) . Pick a small constant 5 > 0, then from u{y/B) = 0, we may 
integrate by parts and compute: 

\\{H b — 5)u\\ 2 L 2 = ({H b —5)u, {H b — 5)u) b = \\H b u\\ 2 L 2 —25(H b u, u) b +5 2 \\u\\ 2 L 2 ■ 

p,b p,B p,B 


We now use the isometry (12.971) and (12.91) . We first obtain from (I2.32|) the spectral 
gap: 

Vu with {u,if B j) B =0, 0 < j < k, 

then 


( H b u , u)b > 


2k+ 4 + 0 


1 


log A? | 


mV , 

p,B 


(B.5) 


and similarly from (1A.ll) : 
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We therefore conclude that for S small enough: 

\\Hbu\H, > \\d z u\\ 2 L2 b + \\{1 + z)u\\ 2 L 2 . (B.7) 

p,B p,B 

We integrate by parts using the general formula 

(d z u, v) B = -x r Bu(Vb)v(\ r B)p(\/ r B) - (u, d z v) B ~ (u, y)e - ( u , zv) B (B.8) 
to compute: 

\\H b u \\ 2 L 2 = \\Au\\ 2 l2 +2{u zz + — ,zd z u) B + \\Au \\ 2 L 2 

P,B u p,B Z U p,B 

= ||Art||^ 2 + (d z (d z u) 2 , z) B + 2 \\d z u \\ 2 L 2 + \\Au\\ 2 l2 

p,B p,B p,B 

= \\Au \\ 2 L 2 - B(d z u) 2 (Vb)p(Vb ) - ((d z u) 2 , l) B - (( d z u ) 2 , l) b 

p,B 

~ \\Au \\ 2 L 2 + 2 \\d z u \\ 2 L 2 + ||Au||i a 

p,B p,B p,B 

= \\Au \\ 2 L 2 - B(d z u) 2 (VB)p(VB). (B.9) 

p,B 

On the other hand, let 


f 1 for \/~B < z < 1 
\ 0 for z > 2 


then integrating by parts: 

(H B u,x(z)\ogz) B = ~{u,H B (x(z)logz)) B + VBd z u(VB)\og(VB) 
and thus: 


B\d z u(VB )\ 2 < 


| logS | 2 


II H b u 


|2 
I L 2 


p,B 


+ M 


"p.Sj 


< — l — 

~ | log S| 2 


\\H b u\ 


L 2 ) 

n p,B 


where we used (IB.71) . Together with (IB.9}) . (IB.71) . claim (IB.21) follows. 


step 2 Proof of (I II. •'ll) . (1B.4D . Define the radially symmetric function 


v(z) = 


( H b u(VB) for 0 < z < VB 


H b u(z ) for 2 > VB 


and note that v € Consider 


(B.10) 


(v,Pj) o 


P; 


w:=v~y ——x— ±n. 

JA'o (Pp P j)o 


Then from (12.111) : 


[ \d* 

Jz> 0 


,2 *1 

w e2 zdz > 


( 2 k + 4) [ | 

Jz> 0 


2 *£ 
b e 2 zdz 


> ( 2 k + 4) / \w\ 2 e 2 zdz. 


On the other hand, from (12.27}) . (IB. ID : 


Pi, ~ 


| logS | 


4>B,k 


z>y/b 


< — i — 
~ llogBI 

p,B 


(B.ll) 
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from which for 0 < j < k: 

(v, p,) 0 = H b u{Vb) 


r ~ z 2 2 f\\ H Bu \\ L 2 ' 

L<vb Pie ' ziz+ {Hbu ■ bj)b+ 0 (nspf" 


/ \\H B u \\ L 2 \ 

= o(siW 5 >i+ nEi 5 f*J. 

where we used the orthogonality {u,if B j) B = 0, 0 < j < k. Therefore 

ll” - to|| h;,b ~ B\Hbu(Vb)\ + 


Injecting this into (IB .111) yields (IB.31) . We now apply (IB. 61) to v and conclude from 

(ED) , (iBTOl) : 


\\zH B u\\ 2 L2 < \\zvf T? _ + j B 2 |F b u(Vs)| 2 < 

p,B 


lL lo 


J Z U \\ L 2 


p,0 


+ ||n|| 2 2 o + B 2 \H b u{Vb){ 


< \\d z H B u\\l 2 +|| J ffRu|| 2 2 +B\H b u{Vb)\ 2 < \\d z H B u\\ 2 L2 +B\H b u{Vb)\ 


" P,B 


- 1 P,B 


P,B 


and (IB.4p is proved. 


□ 


In analogy to Lemma IA.31 we now renormalize Lemma IA.3I by letting e(y) = 
u(y/By) and obtain: 

Lemma B.2 (Coercivity of T~i B ). Let k € N and 0 < B < B*(k ) small enough. If 
£ e nl(v > 1) satisfies 


e(l) = 0, (e,riB,j)b = 0 for 0 < j < k, 


then 


Moreover , 


b ~ \\^ £ \\ 2 b+ b\\d y e\\ 2 B + B 2 \\A.e\\ 2 B 

+ B 2 1|(1 + yfBy)e ||| + B\ log B\ 2 (d y e) 2 (l). 


(B.12) 


%'Ubs || L 2 > 

p,B 


2k+ 4 + 0 


1 


and 


log B | 

1 


£||^e||!-c fc £ 2 |^ Be (l)| 2 (B.13) 


S||yWB£||| < + |^s£(l)| 2 . 


(B.14) 


Appendix C. Non trivial melting initial data 

In this appendix, we show that our set of initial data for melting is non empty 
and contains compactly supported data arbitrarily small in H l . We show the con¬ 
struction for k = 0, an analogous construction holds for k > 1 and is left to the 
reader. 

To see this define a cut-off function x{v) = 1 for y < 1 and x(u) = 0 f° r 2 / > 2, and 
set XB=X{ B )> where 

B 2 = (logbol 


2feo 


By abuse of notation, we denote by bo the initial value of bo in this section. Let 

\\Vb 0 \\b 2 


{XBVbo j hb() ) bo 


a : = 


(C.l) 
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Note that a — 1 = ^ XB)>;b 0 .>7t , 0 )b 0 . furthermore 

(XBVb 0 ,Vb 0 )b 0 


|((1 - XB)Vb 0 ,Vb 0 )b 0 \ < / Vb 0 Pbo ^ / log 2 (y)/) &0 + ||t?6 0 ,iII? 

J y>B J y>B 

r co 


< 


1 f°° 1 

— / log 2 (v^o^)p(g)rfz + 

Oo &o|log6 0 | 2 

log 2 6 0 r , w , 1 r ATT, , 1 

- / p(z)dz + — ^p(z)dz + ———g 

Jv^B ^0 Jv^B fro|log6o| 2 


°0 J^B 

< lo g 2b 0 -6 qB 2 /2 , ]_ -b 0 B 2 /4 , 1 

~ 6 0 6 0 |log6 0 | 2 

< ^°g ^O c -|logfe 0 |/4 J_ e -|logb 0 |/8 _|_ -*■ < 

~ &0 ^0 


6 0 |log6 0 | 2 ~ b 0 |log6 0 | 2 ' 


(C.2) 


where we used (I2.77p and (12.781) . Therefore 

i 


la - 1| < 


bo | logfto | 2 


< 


1 1 
Mlog&ol 2 < 6o|logb 0 | 2 


1 


W'lboWbo ~ (C 1 “ XB)Vb 0 ,Vbo)b 0 ~ Ikollfto ~ l lo g & o| 

where we used (1C.2D and (12.78p . Consider the initial data 

-u(O) = u(0) = b 0 axBVb 0 - 

Note that by (1C.2D we have the bound |a| < 1. Then using (I2.77D . (I2.78D : 

\Vb 0 \ 21 


4 ’ 


(C.3) 


[\dyum 2 < b 2 0 [ | d yVb0 \ 2 +[ 

J Ul<v<2B JE 


< b 2 ^' 2 

< 1. 


1<2/<2B J Bo<y<2Bo 

ll^%ollfeo + ^IKII&O 


B 2 


2 Ji£|M 
^b 0 e 4 


|log6 0 1 + 


bo |log&o| 


K 2b, 


'o 


On the other hand, we have by definition e(0) = — (1 — axB)boVb 0 = — (1 —a)6o%o ~~ 
a(l — XB)boPbo an d hence: 


ll^6 0 e(0)||f <^o I \boh 0 Vb 0 \ 2 pb 0 ydy + \l-a\ 2 b%\ 2 bo \\ri bo \ 

Jy>B 


2 

bo 


+ b o 


L 


\dyVb 0 \ 2 . Pb 0 , ,2 2 


'B<y<2B 

Using (1C.2ft we can estimate the first term: 


2 3 4 + b 0 p b 

y 2 y 4 


Pb 0 ydy. 


b o / \bo\ 0 Pbo\ 2 Pb 0 ydy < b^\ 2 bo < 6qA 

Jy>B 


2 1 < b l 

bo b 0 \logb 0 \ 2 ~ |log£>o| 4 ’ 


For the second term we use (1C.3D and (I2.78p and readily obtain 


|i-«I 2 ^a 2 0 || %0 || 2 0 <^ 


b o 

|log6 0 | f 
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Similarity using the decomposition (12.771) (with k = 0), we have 


'L 


\dyVb 0 ? 


1 


B<y<2B V 


Pb 0 y dy<bl / ^jp bo ydy + b z 0 

JB<y<2B y 


\9yVb 0 , 


B<y<2B 
2 


1 I ^ 

—Pb 0 ydy 


<blB~ A e~ b ° B /2 + b 2 0 B- 2 \\d yVboAUbo 

h 3 

+_< ft o 


b 4 0 


~ llog^ol 4 llogfool 4 ~ |log6o| 4 ’ 
where we used (12.781) in the last line. In a similar fashion 

| logy | 2 


bl 


L 


V bo 

~ipb 0 y 

B<y<2B y 


dy<b 2 0 

JB<y<2B y 

< b 2 B~ A \\ogB\ 2 f 


f 7 

Pb 0 ydy + bl / — '^-p bo ydy 

Jb' — 


>B<y<2B y 
l2 d—4||— 112 


' B<y<2B 


p bo ydy + b 0 B \\y bo ,i\\ bo 


~ llog&ol 2 


M_ e -b 0 B 2 /2 


h 3 

°0 


|log6 0 | e 




~ llog^ol 2 


o_ -|logfeo|/4 


+ 


h 3 

°0 


< 




if is sufficiently small. Finally, 


bl 


/ b 0 V bo Pb 0 ydy <b 0 

B<y<2B JB<y<2B 


|logty| 6 ~ |log6 0 1 6 
(\ogyfp bo ydy+ bl\\p bo ^\\ 2 bQ 


< b^\\ogB\ 2 ^-e boB2/2 + bi 1 

bo &o|log6or 

< fegllog^ope -1 ^+6 q 1 


< 


h 3 

°o 


°|log6 0 | 2 ~ |log6 0 1 2 ’ 

and hence (j3.15)1 is satisfied. Moreover, 

O(0),r?6 0 )& 0 = (-(1 - axB)b 0 y bo ,y bo ) bo = -b 0 \\y bo \\l 0 + ab 0 (xb 0 % 0 ,Vb 0 ) = 0 

by (1C.II) . and therefore the orthogonality condition from (|3.14l) is satisfied. 

Remark C.l. Observe that by our construction, the initial temperature uq is non¬ 
negative in R In this case, the solution u(t, •) remains non-negative by the maximum 
principle. 


Appendix D. Cauchy theory in H 1 x H 2 

Theorem D.l (Well-posedness in H 2 ). Let uo € H 2 (Ll( 0)), Ao > 0, rto(Ao) = 0. 
Then there exists a time T = T(|tyo||# 2 m), Ao) > 0, a constant C > 0, and a 
solution (u, A) to the Stefan problem (11.21) on the time interval [0, T] such that 

u€C([0,T),H 2 (n))nL 2 ([0,T),H 3 (Q)), 

u t €C((0,T),L 2 (n))nL 2 ([0,T},H 1 (n)), 

A€C 1 ([0,T),R), (D.l) 

and the following bounds hold: 

INI <C = C(|tyo||H 2 (a(o)),A 0 )! Hi) > y 
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for some universal polynomial function C of the initial data. Moreover, if T is the 
maximal time of existence of a solution (vj, A) satisfying (ID.II) , then 

either lim ||w(t, OII/DfOdV) = 00 or li m A(t) = 0. 
t->T~ t->T~ 

Remark D.2. The Stefan problem allows for an instant smoothing effect. It is 
well-known that the solution u becomes infinitely smooth on ( 0 , T) in both the 
time- and the space variable, see for instance [26;, [39] [25] . 

The proof of Theorem ID. II is presented at the very end of this section, as a simple 
consequence of Theorem ID.51 

We start by pulling-back the problem (11.21) onto the fixed domain If := {y € 
R 2 , |y| > 1}- We denote the points in D by bold y, while the radial coordinate |y| 
is denoted by y. We define the pull-back temperature function w : 12 —)• R by 

w{t,y) = u{t,\(t)y) 

A simple application of the chain rule gives the following system of equations for 
w : 


A 

A" 


— —7 Aw = 0 in 12 ; 

A^ 

(D.2a) 

w y (t, 1) = -A(f)A(f); 

(D.2b) 

w(t, 1 ) = 0 , 

(D.2c) 

II 

S 

O 

V 

II 

V 
0 

(D.2d) 


Lemma D.3 (Energy identities). Assume that (w, A) is a smooth solution to the 
Stefan problem (ID.21) on some interval [0, T\. Assume that A(0) > 0, tco |, y _ 1 = 0, 
and that w(t,-) € H 2 (Q) for t € [0, T]. Then on the interval [0,T] the following 
energy identities hold: 


9^IMIi 2 (Q) + ^ 2 ll v ' l ^’llz, 2 (o) - ~IMI! 2 (q)> ( d - 3 ) 

^llVHIi^n) + ^\\ Aw \\l^m = ^A 3 , (D.4) 

\j t \\ Aw W 2 L^a) + YJt ( A I A 0 + 5(2 l! VA ^ll^(n) 

9 A 

= -j\\ Aw \\h(n) + -A 5 A 3 . (D.5) 

Proof. Multiply (ID.2a.l) by w and integrate over Q. We obtain 

0 = hit i “’ 2 rfy ~ 1 / n A "™ dy +X 11 ' Vw|| l j <"> 

1 d „ 112 A f°° , on , 1 1,2 

= 2 dt^ L2 W “ 2A y dy ^ w h d y+^\\^ w \\l^i) 

= 2 ^IHli 2 (n) + ^IMli 2 (n) + ^2 ll Vu; ll i 2 (n)> 
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which is precisely (ID. 31) . Multiply (1 1).2a I) by —Aw and integrate by parts. We 
obtain 

0 = J Vw ■ Vw t dy - 2Tr(d n ww t )\ y=1 + jJ AwAwdy + -^ || Arc]|| 2(n) 


1 d 


A 


1 


2dt W s/w Wh(n) +‘lTT{w y w t )\ y=1 + 27r-y Awd y Awdy + ^ || Aw|| 2 2(n) 

^H Vu; lli2(a) - 7r ^( Au; ) 2 | y= i + y W Aw \\lw 

1 d ,,2 


1 


w. 


y\ y =i 


~ 2dtW^ W ^ L2 ( n) ~ + W ll Au; llL 2 (f2)> 

where we used the fact that wt(t, 1 ) = 0, Arc = ^d y Aw, and Arc|^_ 1 = 

—AA. Note also that d n w\ y= \ = —9 y rc| y= i.This proves (ID.41) . To prove (ID.5|) we 
apply V to (ID. 2D . multiply by — VA w and integrate-by-parts. We obtain 

0 = — [ \7w t ■ VAwdy + ^ f VArc • VAwdy + ||VAic||^ 2 (q) 


in 


\ d f \ f 

--j~ t j (Aw) 2 dy + 2n(dyWtAw)\ y=1 — — J AAw-Awdy 

- 2irj(d y AwAw)\ y=1 + ^ ||VAw||| 2(n) . 


(D.6) 


From (lD.2bl) it follows that 


dyw t \ y=1 = ~d t { AA) 


Restricting (ID.2al) to y = 1 we conclude that —jw y \ y _ 1 — p'(Arc)| y _ 1 = 0. Us¬ 
ing (ID.2bp this implies that 

(Aw)| y=1 = A 2 A 2 . (D.7) 

The previous two boundary identities imply that 


2ir(d y w t Aw)\ =-^-^(aA) . 


'y=t 3 dt V"7 ' ^ D ' 8 ^ 

Note that d y Aw = w y + yw yy and therefore, when restricted to y = 1 we conclude 
that (c^Aw)] = (Aw)| . From (ID.71) we infer that 


(d y Aw)\ = A 2 A 2 . 


Therefore 


A 


— 27t— (d y AwAw)\ r =1 = —27tA 3 A 5 

a y 


(D.9) 


It remains to evaluate the term — j/q A Aw ■ Awdy. A direct calculation yields 

AAw = AAw + 2Aw. 

Therefore 

\ r 2A A 

— — J AAw-Awdy = — —\\Aw\\ 2 L 2^—2 tt— J d y AwAwdy 


= -y l|Am|| 2 2 ( Q ) + TTj{Aw) 2 \ y=r 
Plugging (ID. 81) . (ID.9p . and (ID.101) into (ID.61) . we obtain (ID.5p . 


(D.10) 

□ 
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Let us define the energy-like quantities 

E (t) = o SupJ^||w(s, Olli^fi) + ^||Vw(s, -)llia (n) + ^]| Aw(s, 01112(0)}, (D-ll) 


D ( t ) ~ 'x(fj 2 ^ w ^ t, Ollx 2 (Q) + 3^2 0)11x2(0) + ||VA«;(i,Ollx 2 (o). 

(D.12) 


Lemma D.4 (A priori estimate). Assume that (w,\) is a smooth solution to (ID.2D 
on some interval [0, T*]. Assume that Ao > 0, ioo| j = 0, and that w(t , •) € H 2 fA) 
fort £ [0, T*]. Then there exists aT = T(E{ 0),Ao) > 0, T < T*, such that for any 
t £ [0, T] the following a priori bounds hold: 

Eft) < 4£(0), (D.13) 

A(t) > y- (D.14) 

Proof. From Lemma lD.31 we infer that 

E (t) + J D(s)ds <E( 0) + —AqA(O) 3 ——A(i) 3 A(i) 3 + J 


IAQQI 

A(s) 


ds Eft) 


+ vr 


/' 

Jo 


A 3 (s)|A(s)| b + A(s)A(s) 3 ) ds 


(D.15) 


Note that 

Therefore 


w y\ y =i ~ l|Aw|lx 2 (n) + l|V^||| 2( o). 


\X(t)\ 2 <^E(tf C>1. 


(D.16) 


Using (ID. 16|) we obtain 

rt |A(s)| 


E ®J dS + 7T J 0 ( A3 ( S )|A( S )| 5 + A(s)|A(s)| 3 ) ds 

< C't (E 2 {t) + E 3 / 2 (t) + E 5 / 2 (t)) sup 

V ' o<s<t AUs) 


(D.17) 


To bound the error term — 4^A(t) 3 A(f) 3 we need a more refined estimate than (ID. 16|) 
due to the absence of the integral-in-time. Note that by the trace inequality and 
the interpolation between fractional Sobolev spaces, we have 

K( 1)1 < 11 Vu; 11^1/2(0) < HV^||i / 2 2 (n) ||Vu;||^ / 1 2 (n) . 

Therefore, upon using the Young inequality 

|AA| 3 < ||V«;||^ 2 (n) ||Vti;£ / 1 2 (n) < S\\Vw\\ 2 Hl{n) + C s \\Vw\\ 6 L2{n) 

<6E + C s \\Vw\\ 6 L2{n y (D. 18 ) 

Integrating (ID.4|1 over [0, t]. we have 

sup ||VH|!2 (n) < l|V^o|li 2 ( o) +vr [ A(s)|A(s)| 3 ds 
0<s<t Jo 

< E 0 + C*tE 3/2 (t) sup -- 27 -r. 

o <8<t A A [s) 
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Therefore, we obtain from (ID. 181) that 

|AA| 3 < C 0 + 6E + Ctp(E)q( sup 1 ) 

o <s<t A z (s) 

where p and q are increasing polynomial functions of their arguments. Plugging this 
bound back into (ID.15p . using (ID.171) . and the definition of E[t) we conclude that 

E(t) + [ D(s) ds < C 0 + Ctp(E(t))q( sup 1 ), (D.19) 

Jo o <s<t * (s) 

where Co = Cq{E{ 0), Ao). Since A (t) = A(0) + J Q £ A(s) ds, it follows that 

Ct 

A (t) > Ao — t sup |A(s)| > Ao-— E l ^ 2 {t). 

o <s<t A 


Let 

T' = sup{f > 0| E(t) < 4E( 0), A (t) > Ao/2}. 

By the continuity of E(-) and A(-) it follows that T' > 0. On [0,T X ] we therefore 
have 

E{t) + f D(s) ds<Co + Ctp(E(t))q( - ^ --) (D.20) 

Jo Ao -4^§f£(0)V2 

By a standard continuity argument, there exists a sufficiently small T = T(E( 0), Ao), 

T < — _'*° , such that 

— loVCEio) 1 / 2 

E(t) < 2(7 0 , t G [0, T], 

By the choice of T, we also have the bound A (t) > ^Ao > ^Ao for t € [0, T] and this 
concludes the proof of the lemma. 

□ 

Theorem D.5 (Local well-posedness). Let wq G H 2 (Ll), Ao > 0, and u;o| ( _ 1 = 0. 
Then there exists a time T = T(||u>o||H 2 (r 2 )> A(0)) > 0 and a solution (w,X) to the 
Stefan problem (El) on the time interval [0, T] such that 

w € c([o, t],h 2 (Q)) n l 2 ([o,t], h 3 (Q)), 
w t eC((0,T},L 2 (n))n L 2 ([0, T], H 1 (O)), 

A G C' 1 ([0,T],R), (D.21) 

satisfying the energy estimate 


E(t)<Co = C 0 (E(0)), t € [0, T] 


and the lower bound 


X(t) > 


A(0) 

2 


t € [0,T], 


where the energy E(-) is defined by (ID.111) . Moreover, if T is the maximal time of 
existence of a solution (w, A) satisfying (ID. 211) . then 


either lim llWi, •)||// 2 (a'i = °o or lim X(t) = 0. 
t-s-T- t^T~ 


Proof. The proof of existence follows a standard iteration argument for the sequence 
of approximations (X n (t), w n (t)), n € N. For a given A n (-) we define w n +i by solving 

- '^ !i Aw n+ i - -^-Aw n+1 = 0 in Ll ; 

An X n 

w n+ i(t, 1) = 0. 


d t w n+ 1 


(D.22a) 
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We then update A n +i(-) by solving 

9yVJn+i{t,l) — dt A n -|_i (i)A n (t). 

Estimates analogous to the a priori estimates of Lemma ID. 41 can be used to ob¬ 
tain uniform bounds on E(w n , X n ) + D(w n , X n ) ds, where E and D are defined 
by (ID.Ill) and (ID.121) respectively. Note that we can also get uniform bounds 
on ||<9tWn||i,°°([0,T],L 2 (n)) + ||^n||L2([0,T],i?i(n)) as the latter norms are controlled 
by E(w n ,X n ) + Jq D(w n , X n ) ds from (ID.22ap . Upon passing to the limit, we ob¬ 
tain a solution to which the energy estimate of Lemma ID.41 applies. The proof of 
uniqueness is standard. The breakdown criterion is a simple consequence of (ID. 111) , 
and (ID. 161) . □ 

Proof of Theorem ID.ll Let w be the solution to (ID. 21) as given by Theorem ID. 51 
It is easy to check that \\w(t, •)IIh 2 (o) 1 $ t £ [0,7~). Theorem ID.ll now follows 
from the change of variables u(t,r ) = w(t, jfo) and Theorem ID. 51 
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